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Result*  frcss  tfee  Hmse  H.Ung-BiawtlMi  fii*h“£xs?i©si¥«  SSsHistiea  YecJsaliSiir 
CLDSSI5TJ  exrer i&est  are  presented .  fM®  eaperlgest  sia&Iatess  the  air-MASt- 
instated  ground  aotiem  free  &  large-yield  nuclear  ©eape*.  over  &  plan  area 
94  feet  fey  150  feet,  k  specially  designed  aetrix  of  prUsaeord  «aa  utilised 
to  produce  the  desired  espies  lets  is  &  eoaf  ined  volua©  ©f  air.  ever-fjordas* 
ar  ,mss  of  ss&ierial*  was  piacoi  ©ver  the  expi^stee  t©  pwi4e  a  rg&ciisi®; 

T&tt.ti  to  shape  tiae  ressdtist  wave  pals©  and  to  provide  the  repaired  is®g 
durations.  fhe  &xper  least  produced  as  air  pressure  pula©  afeich  had.  a  ptmk 
©verjsressa r#  of  312  psi,.  e  tige  to  fiee-helf  peak  pressure  of*  IS. 2  asaec »  asl 
a  total  duration  ©f  1T0  aaec.  fhe  stsefe  frsnt  traveled  at  a®  *srgr»«e  vtlmitj 
©f  $,X2$  ft/sec.  S®  -the  feasts  of  theee  rossdLis,  r«^Mfesa^ati«^  ©re  sale  for 
future  riastiaii&s  experiments .  the  assured  saris  #ra®-fi«M  asstiosss  ana 
stresses  are  e-oapered  vith  theoretical  saleulaticns,  and  the  validity  of  the 
technique  for  slisuiaUng  the  desired  nuclear  envir^rateet  is  estahLiafeei . 
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The  P&ase  II,  Umg-&*u  ratio®  Migh-SEpiosive  Slsssiialios  Technique  C LAMEST ] 
shot,  was  fired  at  1315  boar#,  i©T,  on  1$  ©ecejaber  19^»,  at  Kirtland  Air  Force 
Base,  Mew  Mexico. 


The  primary  objective  of  this  test  wee  id  eiaulat*  the  air-hlast-iadsioed 
grrnisd  actions  fress  a  iarge~yiei«S  nuclear  surface  explosion.  T&e  specific 
reipireesest  was  to  siMulatsr  the  eavircsssasEt  at  tsc  3aO»psi  overpressure  r&sg® 
from  a  large-yield  mieleajr  explosion.  Secondary  objectives  wear#  to  observe 
the  response  of  a  scaled  structural  ss&del  aad  to  proof  teat  a  hardened  cable 
complex  in  the  stsuiaied  environment . 


t 
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2.  Background 


A  great  deal  of  interest  has  been  recently  generated  in  the  development 
of  various  techniques  to  simulate  nuclear  weapon  effects,  tteeh  of  the  impetus 
for  this  increased  interest  resulted  fro®  the  saoratorius  on  atsscesf&eric  testing 
which  was  imposed  fey  the  Sue 1 ear  Test  Ban  Treaty.  Air  Force  requirements  exist 
to  conduct  nuclear  tests  or.  actual  hardened  systems  as  well  as  to  conduct 
basic  research  ©s  miclu&r  weapon  effects.  Since  atmospheric  testing  is  so 
longer  possible,  these  requirements  must  now  fee  act  through  simulation  prograss®. 
To  meet  a  portion  of  these  requirements,  the  Air  Force  Weapons  Laboratory  began 
a  two-phase  progrsie  in  February  19^^  to  develop  a  technique  to  simulate  the 
air-biast-induced  ground  jactic-ss  fro®  &  large-yield  midear  explosion.  Phase  I 
(reference  1}  of  this  program  consisted  of  a  thorough  study  of  various  simula¬ 
tion  tec  snivel;  and  the  actual  field  testing  cf  the  two  aost  promising 
explosive  techniques.  The  techniques  tested  were  {1}  a,  deton able  gaseous 
mixture  and  (?)  a  matrix  of  prisaccrd  to  produce  the  desired  explosion  ic  a 
confined  volume.  An  overburden,  or  mass  of  material,  was  placed  over  the  ex¬ 
plosive  tc  provide  &  reaction  force  which  would  snap*  the  resulting  wave  pulse 
tc  provide  the  requires  long  durations. 
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Tm  ajspl  icatisG  of  &r<  overburden  tc  increase  explesiee  durations  was 
first  attempted  'ey  the  Stanford  8« search  Institute  (reference  2} .  They  fired 
a  saafeer  of  simulation  snots  on  a  small  scale,  5-1/2  feet  x  11-1/?  feet  is 
piss*  utilising,  a  dslusashlt  ^.seoas  sixtare.  t3ss.ee  prelj&isary  results  indi¬ 
cated  that  the  overburden  concept  for  increasing  durations  was  a  practical 
sc  naso.e  if  the  ex  per  iment  was  conducted  on  a  large  enough  scale,  however,  the 
gaseous  mixtures  which  they  utilised  did  not  simulate  the  traveling  wave 
characteristics  of  a  shock  wave  in  air,  i.e.,  the  correct  pressure/f rental 
velocity  relationship#. 

the  pr  Leaner!  matrix  concept,  also  utilising  the  overburden  to  increase 
durations,  was  developed  under  contract  (reference  3)  for  the  Air  Force 
Weapons  Laboratory.  Two  shots  were  conducted  cs  an  area  15  feet  x  30  feet  in 
plan,  these  efforts  were  not  successful  in  producing  the  desired  environment 
of  a  303-psi  overpressure  wav#  from  a  large-yield  nuclear  weapon. 

Seven  shots  were  fired  under  the  Phase  I  study  on  an  area  20  feet  x 
**0  feet  in  plan,  these  tests  farther  developed  the  two  explosive  techniques 
and  refined  the  overburden  concept  for  increasing  durations.  Either  technique 
appeared  to  be  capable  of  Simla  ting,  the  desired  environment,  however,  the 
priaaconi  aatrix,  officially  designated  the  Long-Duration  Eigh-Ex plosive 
Simulation  Technique  {LSEEST),  was  selected  as  the  -acre  promising  for  appli¬ 
cation  on  larger  scales.  It  was  easier  to  field,  safer,  and  aware  economical 
than  the  de tenable  gas  technique. 

The  Phase  II  test  was  a  direct  follow-on  of  the  Phase  I  effort.  This 
test  was  conducted  to  demonstrate  the  feasibility  of  simulating  air-blast- 
induced  ground  actions  from  a  large-yield  nuclear  weapon  utilizing  the  LDEEST. 
The  test  was  conducted  on  a  plan  area  96  feet  x  150  feet,  which  was  much 
larger  than  any  experiment  previously  conducted.  Conventional  construction 
techniques  and  practices  were  utilized  to  demonstrate  the  feasibility  of  con¬ 
ducting  even  larger  experiments. 

3*  Suclear  Blast  Environment 

This  experiment  was  planned  to  simulate  the  nuclear  blast  environment 
at  the  300-psi  peak  overpressure  level  from  a  1-kfT  surface  explosion.  This 
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peak  over  presort*  will  ©gear  at  a  grssaad  range  &t  2,55S  feet  from.  the  detonation 
point  and  will  have  s.c  associated  shock  front  velocity  of  fe,BOO  foot  per  second 
C reference  A),  The  fi^erpreiis^F*:  decays  rapidly  fro®  the  peak  value  at  the  shock 
frost  and  reaches  a  value  of  15-0  p§i  ih  37  sa-ec.  lowsver,  b&e  total  positive 
P&ase  duration  will  be  appro* imately  1  second . 

figure  X~i  illustrates  the  pressure-tise  history  of  &  3CK>~psi  fssak  over¬ 
pressure  shock,  wave  from  a  I -iff  surface  burst.  The  isyul  se  contained  ie  tfels 
pulse  eaa  be  computed  by  taking  the  time  integral  of  the  pressure  over  the 
entire  positive  phase  duration.  In  this  instance*  it  is  32  psi~s*c. 

Mfeen  the  history  of  the  stock  wave  is  studied  l  figure  1-2)  it  can  be  seen 
that  the  peak  overpressure  and  the  shock  front  velocity  Mtfc  decay  •aooetonicaily 
with  ground  range,  This  is  caused  by  the  expansion  taking  place  as  the  air 
blast  loads  increasingly  larger  disks  of  material. 

The  curvature  of  the  shock  front  also  decreases  with  ground  range.  At  a 
distance  of  2,000  feet  from  the  detonation  point,  the  curvature  of  the  shock 
front  Is  considered  to  be  insignificant.  Therefore,  this  simulation  technique 
utilises  a  plane  wave  front  tc  simulate  the  actual  case. 

The  air-blast  parameters  that  were  considered  in  the  development  of  the 
environment  to  be  simulated  are 

Cl)  Peak  overpressure 
(2)  Shock  front  velocity 

v3 )  'Overpressure  duration,  pulse  shape  and  impulse 
(b)  Peak  overpressure  and  shock  front  velocity  decay 

The  effects  of  these  parameters  on  the  air-blast-induced  ground  motions 
are  discussed  in  the  next  section. 

b.  Theory 

a.  Air-Blast -Induced  Ground  Motions 

When  large  surface  areas  are  loaded  nearly  uniformly,  such  as 
occurs  at  moderate  ranges  from  a  large-yield  nuclear  detonation,  one- 
dimensional  wave  theory  can  be  used  as  an  approximation  -  T«i  one— dimensional 
wave  theory,  only  vertical  motions  occur  and  simple  relationships  can  be  es¬ 
tablished  among  stress,  displacement,  strain,  velocity,  and  acceleration. 
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Although  this  situation  is  highly  idealized,  the  theory  previa**  tm  w* Aar*- 
stand  lac  of  the  basic  relationships  between  the  air-felast  parameters  and  the 
stresses  and  actions  induced  in  the  soil.  These  results  can  then  be  corrected 
to  tale  account  &£  the  other  factors  ehi«h  influence  the  free-ixeid  behavior. 


Fiji ire  1-3  shews  a  uniform  elastic  medium*  extending  to  an  infinite 
depth,  which  is  loaded  uniformly  at  it*  free  surface  by  a  tine  dependent 
pressure,  Kt).  A  unit  cylinder  of  soil,  extending  fro*  the  surface  tc  as 
infinite  depth,  can  be  examined  and  sill  be  representative  of  all  the  sur¬ 
rounding  material.  The  strains  in  all  directions  except  depth,  t,  are  equal 
to  zero. 

For  equlllbricsi  the  summation  of  forces  in  the  x  direction  must 
he  equal  to  sere,  and 


(pdxl 


ds  »  0 


which  simplifies  to 


o 


3^U 

it7 


where 

c  *  mass  density  of  the  medium 

s  «  undisturbed  location  of  a  given  particle 

u  *  displacement  of  a  given  particle  in  the 
s  direction,  a  function  of  *  and  t 

t  *  time  variable 

ax  *  normal  stress  in  z  direction,  tension 
considered  positive 

Define  strain  as 


a  •  M< 

x  z 


and  stress  as 
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where 


*  coostraises  msdul us 
£  *  £ssag?s  modulus 

v  *  Poisson's  ratio 


Then  making  appropriate  safe stitmticns  into  the  wave  equation. 


«  EHIm 
"s? 


Sefine 


€2  *  K/n 


fber.  the  wave  equation  reduces  tc 

i2^  _  «2 
at?  *  c  li? 

The  solution  to  this  equation  is  of  the  form 

u  *  f{t  -  x/C}  ♦•  g{t  ♦  z/C) 

where 

fCt  -  z/C)  describes  a  forward  mowing  wav* 
g{t  ♦  x/C)  describes  a  backward  mowing  ware 

Since  the  medium  is  considered  to  hare  an  infinite  depth,  only  the  forward 
•owing  ware  solution  is  required. 

The  ware  shape  is  unaltered  as  it  propagates  through  the  medium. 
Figure  X-i  shows  a  pressure  pulse  propagating  through  the  medium  with  the  time 
variable,  (t  -  x/C).  z/C  is  the  time  it  takes  the  ware  to  propagate  from  the 
surface  to  the  depth  being  studied. 

Stress  has  been  previously  defined  as 
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with  appropriate  substitutions 

o  *  C*ct  *  C*p  *  H>Cf*{t  -  */C) 

£  X  #1 

vticre 


u  ■  fit  -  */c) 

*  =g-  rtt  -  */c) 


strain  has  been  previously  defined 


►  .  3a 

*  as 


which  free  above  is 


tz  *  f*{t  -  */C> 
Particle  velocity  Is  defined  as 

u  *  ft  *  *'  ^  ~  i//c' 

At  tee  surface  (*  *  0)  the  stress  is  equal  to 

az  *  -P(t) 


Then  froa  above 


at  «  -P{t)  «  -oCf‘(t  -  x/C) 
and 

f*{t  -  x/C)  • 

At  a  depth ,  t  f‘(t  -  x/C )  •  P-^- 
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Making  the  appropriate  substitution* 

o  «  -?(t  -  */C) 

A 

£  »  -?(t  -  */C}/pC2 

1 

u  »  ?{t  -  */C)/tfC 

Acceleration  is  defined  as  the  parti  si  derivative  of  the  parti  cle 
velocity  with  respect  to  tise 

«  *  «  -  ^  Him  -  */C)) 


Thus  ,  the  peak  acceleration  is  dependent  on  the  shape  of  the  over* 
pressure-tiae  curve. 

The  absolute  di  splaceaent  of  a  point  *  z ,  at  a  specific  tiae,  t  ,  can 
he  obtained  by  art  integration  of  the  particle  velocity  at  the  point*  z 


uiz,  t)  •  £§  J 


T{t  -  z/C)dt 


and 


a(z,  t)  « 


where 


I 


P(t  -  z/C)dt 


Tais  integral  is  equal  to  the  area  under  the  overpreaaure-tiaa  curve 
P(t  -  z/C;  between  the  tiars  o  to  t.  See  figure  I-h.  This  integration  can 
he  conveniently  perfomed  graphically  on  the  overpresaure-tiae  curve. 

Table  1*1  suaaarizes  the  one-diaenaional  wave  propagation  quantities. 
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Table  1-1 

msmm  or  tm-mmssszomL  mm  raorASATios  wms 


quantity 

Equation 

Sire** 

*  -P(t  -x;€) 

Strain 

c  *  -F{t  -  s/C}/p€5 

Velocity 

u  *  P(t  -  i/C) /oC 

Acceleration 

*  *  £c  at  {Ht  -  I/C)J 

Displacement 

u  *  I/pC 

The  one-rU  —tie  too  el  wave  theory  has  shorn  that  the  stress  in  the  soil 
propagates  with  a  ware  shape  unaltered  from  the  air-blast  input .  In  reality, 
three  major  changes  will  take  place  as  the  stress  ware  propagates  with  depth. 
The  rise  time  will  increase,  the  impulse  duration  will  be  spread  out,  and  the 
peak  stress  will  decrease.  nonlinear! ties  in  the  stress-strain  curses  of 
earthen  Materials  are  the  priaary  cause  of  the  first  two  factors.  The  varia¬ 
tion  of  stress  intensity  with  depth  is  caused  by  three  separate  phenomena. 
First,  aost  earthen  aste  rials  have  highly  none  last  ic  stress-strain  curves 
which  result  in  energy  absorption.  Secondly,  when  finite  areas  are  loaded 
by  transient  Inputs,  rarefaction  wares  propagate  inward  froa  the  boundaries 
(see  figure  1-5)-.  Thirdly,  spatial  attenuation,  or  the  spreading  out  of  the 
load,  and  reinforcing  of  the  load  by  areas  that  have  been  loaded  by  higher 
pressures  play  conflicting  roles.  As  previously  described,  the  peak  over¬ 
pressure  and  shock  velocity  of  the  air-blast  input  are  both  decaying  with 
ground  range  (see  figure  1-2).  If  this  air-blast  Is  used  as  the  loading 
input  for  an  elastic  half  space,  ray  theory  can  be  used  to  exaelne  the  problem 
(see  figure  '  ~6) .  Each  peak  overpressure  le^el  has  a  r*y  path  along  which 
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is.  propagates  ,  and  spatial  41  verges®*  causes  it  to  attenuate  as  it  propagates, 
if  &  vertical  section  is  examined ,  it  can  fee  shown  that  energy  fron  higher 
overpressure  regions  will  fee  feesiisg  late  and  locally  reioferei&g  the  stress 
level  at  any  range  of  interest,  For  negato®  yields,  the  sir-blast 

parameters  are  changing  at  the  appropriate  rate  to  balance  the  effect  at 
the  spatial  divergence  amd  »o  attesnatifim  #j@  to  this  ptfeesaMsefioss  nhonid  ©ee^r. 

s fetes  the  teat  area  is  small  C *'500  feet)*  in  comparts-em  with  the  area 
loaded  ey  a  i-MT  weapon  *  the  air-blast  parameters-  are  nearly  mi  fares  (see 
figssn?  1-2).  In  sd&ities,  ?<feea  depths  ©f  a  fsw  hundred  feet  or  less  are  of 
ict*r»st,  the  ruling  effect  is  minimi.  Because  of  these  factors ,  tMe  decay 
of  peak  overpress  are  and  shock  front  velocity  bairn  not  been  farther  considered 
as  tit  important  parameter  to  the  simulation  technique . 

The  one  -  dimens  i  onal  nave  theory  has  also  shown  that  the  strata  and 
particle  velocity  are  affected  by  the  magnitude  and  wave  shape  of  the  air- 
blsfrt  input  ana  the  properties  of  the  medium.  These  free  field  quantities 
are  also  affected  by  the  sale  parameters  which  cause  the  attenuation  of 
stress  with  depth.  Acceleration  has  been  previously  shown  to  be  extremely 
sensitive  to  the  shape  of  the  air-blast  input,  particularly  the  rise  time, 
which  increases  with  depth.  Displacement  has  been  shown  to  be  dependent  an 
the  integral  of  the  stress-time  history  at  the  depth  of  interest .  fe 
achieve  significant  displacements ,  the  impulse  must  grossly  match  that  of  the 
nuclear  explosion  and  the  area  loaded  mast  be  significantly  large,  in 
relationship  to  the  depth  of  interest,  to-  minimite  the  effect  of  rarefaction 
waves. 

Only  vertical  components  have  been  discussed  thus  far.  A  traveling 
wave  must  be  considered  tc  study  the  relationship  between  vertical  and 
horizontal  stresses  and  act  ions.  Figure  1-T  Illustrates  a  typical  air-blast 
input.  Pit)  ,  traveling  over  the  surface  of  the  earth  with  a  shock  velocity,  U. 
The  velocity  of  propagation  of  the  stress  wave  into  the  earth,  C,  is  deter¬ 
mine  a  by  the  intensity  of  the  air-blast  input  and  the  properties  of  the  earth 
material.  It  is  given  by 


m  mm"  ■ 

c  «  4t7fT 
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when 

If  *  coos  trained  modulus  at  the  stress  level  of  interest 
p  •  mss  density 

For  the  sake  of  simplicity  only  the  si  shiest  input  is  the  superaelanic 
region  will  he  discussed,  Shis  is  s  good  assusption  for  facilities  which  are 
constructed  t *  soil  or  soft  rod  to  withstand  overpressures  of  several 
hundred  p*£ .  In  this  region,  the  stress  wave  trails  behind  the  air-blast 
shod  front  and  enters  the  ground  with  an  angle,  a,  which  can  he  determined 
by 

sin  a  *  C/U 

As  indicated  in  figure  I-?,  the  initial  stresses  and  actions  propagate 
in  a  direction  perpend! cular  to  the  stress  wave  front.  However,  behind  the 
wave  front  these  directions  change.  At  very  late  time*  the  problem  can  be 
considered  static  and  the  motions  and  stresses  must  now  be  vertical,  the 
direction  of  the  motions  and  stresses  thus  change  with  tine,  having  an 
initial  direction  perpendicular  to  the  compression  wave  front  and  a  final 
direction  on  a  vertical  line.  As  an  approximation,  the  components  of  the 
motions  and  stresses  can  be  computed  by  assuming  that  the  entire  stress  vase 
propagates  in  a  direction  perpendicular  to  the  compression  wave  front.  For 
this  assumption  the  horizontal  components  and  vertical  components  are 

Horizontal  component  *  value  perpendicular  to  front  x  sin  a 

Vertical  component  *  value  perpendicular  to  front  x  cos  a 

Thus ,  the  value  of  the  horizontal  and  vertical  components  of  both 
the  initial  stresses  and  motions  induced  in  the  soil  are  a  function  of  the 
material  properties  as  well  as  the  shock  front  velocity,  to  achieve  appro- 
pri ate  horizontal  and  vertical  components,  the  shock  front  velocity  produced 
by  the  simulation  technique  must  closely  duplicate  that  of  the  nuclear 
environment  at  the  overpressure  level  of  interest. 

The  previous  discussion  of  sir-blast- induced  ground  motion  has 
shown  that  the  peak  overpressure,  shock  front  velocity,  overpressure  duration 
pulse  shape,  and  impulse  are  significant  parameters  la  defining  the  environ- 
neat  to  be  simulated.  The  Phase  IX  test  attempted  to  generate  s  peak  over¬ 
pressure  of  300  pel  and  a  shock  front  velocity  of  i  ,800  ft/sec.  The  duration 
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to  cue -half  fc  pressure  was  m  %  criteria®  f@r  4®  fig  lag  the  e©r  feet  ness 

of  the  overpressure  decay  behind  th«  sho^t  front.  Is  addition,  lb*?  gross 
saiehisg,  of  the  ispiiise  ^fts  attempted. 

b.  Priaacstrd  Hatrix 


Figure  1-6  Illustrates  the  general  eo&cept.  of  ifeo  priisaofflrd  n&lrix. 

A  cootiauous  strard  of  prissheord  has  seen  laced  1st®  a  matrix  sbieh  appml- 

sates  the  usifera  properties  of  &  solid  sheet  explosive,  this  ^pprcxisaiioti 


sine  sot  currently  eosawrsially  aval  Xa&ie .  Prictacc-rd ,  an  the  other  hand,  has 
the  appropriate  detonation  properties  when  it  is  placed  in  a  specially  deal gn*4 
sa&rix.  la  addition,  it  is  csMacrel ally  available,  fairly  eeattestical,  and 
safe  to  work  with. 

there  are  two  peraweters  which  say  he  adjusted  to  control  the  character¬ 
istics  of  the  shock,  wave  produced  fey  this  technique.  These  parameter®  are  the 
load  density  md  the  prinacerd  angle.  To  understand  the  effect  of  these 
two  pansaseters ,  the  process  that  produces  the  shoes  ware  stust  be  considered. 

The  wrap  angle  determines  the  rale  at  which  the  ccaeusticci  product*  are 
forked  along  the  length  of  the  facility.  These  coatustion  products  act  like 
a  piston  which  loads  a  cylinder  of  air  in  front  of  the  detonation.  This 
process  foras  a  shade  ware  ia  the  air  which  sores  out  ahead  of  the  platan* 
and  the  particle  velocity  of  the  air  behind  the  shock  wave  is  equal  to  the 
velocity  of  the  piston.  The  p  nobles  is  further  coop! leaked  by  the  intense* 
short -duration  shocks  that  are  produced  by  the  detonation  of  the  individual 
strands  of  pristacord  and  the  reflection  of  these  shocks  with  one  another  and 
the  boundaries  of  the  cavity.  These  secondary  shock*  wove  forward  at  a 
faster  velocity  and  overtake  the  aaio  shock,  the  net  effect  ia  a  rt  in  force - 
sect,  of  the  sain  shoe*  with  a  corresponding  increase  in  pressure  and  shock 
velocity.  The  interaction  of  shocks  creates  a  natural  pressure  decay 
i Mediately  behind  the  shock  front.  If  the  cavity  were  perfectly  insulated 
and  had  rigid  boundaries ,  the  gases  inside  would  eventually  cone  to  sons 
equilibrium  pressure .  however,  the  overburden  begins  to  aaove  upward  at  early 
tines,  which  increases  the  cavity  voluse  and  causes  a  corresponding  decrease 
in  pressure. 
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Cl!  Peak  Shook  Pseasu se 

The  t&eoreti  cal  prediction  of  the  peak  shock  pressure  proved  to 
be  intractable.  Therefore  ♦  m  experimental  approach  was  used  in  Phase  I 
C reference  1)  to  determine  the  peak  pressure,  This  approach  was  guided  by  a 
theoretical  the rnccheai cal  analysis  where  iseatropic,  equilibrium  pressures 
were  computed  for  a  fixed  volume  Preference  3).  these  calculations  considered 
the  ooafomtioo  of  the  primaccrd  and  its  case  and  the  dis association  of  the 
©cshuktic®  presets.  Figure  1-9  pres-eabs  a  logarithmic  least -square  fit  of 
the  four  experimental  data  points  which  were  available  froa  previous  expert - 
aec ts-  fbe  peak  shock  overpressure  fete  bees  plotted  m  a  fdsetica  of  load 
density.  Load  density  is  defined  as  the  total  weight  of  explosive,  in  pounds, 
divided  by  the  initial  confined  volume  of  air*  in  cubic  fact.  According  to 
figure  1*9,  a  load  density  of  0.076  lbs /ft  is  required  to  produce  a  peak 
shock  overpressure  of  300  psi . 

(2)  Shock  Trent  'Velocity 

Figure  I-c  also  illustrates  a  plane  shock  front  propagating  in 
the  desired  direction  of  travel,  the  primacord  has  a  detonation  velocity  far 
in  excess  of  the  fc.ft 00  ft/sec  shook  front  which  is  to  be  simulated.  therefore, 
the  priaaccrd  matrix  must  be  designed  to  have  a  wrap  angle  which  produces  & 
shock  front  having  the  appropriate  velocity.  Staple  trigonometric  theory 
suggests  that  the  wrap  angle,  ®,  should  be  defined  by  the  following  foraula: 

i*  f  *  shock  front  velocity 
5  L  detonation  front  velocity 

xxperi  mental  data  from  Phase  I  indicate  that  this  simple  theory 
does  not  explain  the  observed  results,  the  problem  is  coaqplicated  by  all  of 
the  factors  discussed  in  section  X.fc.b.  therefore,  the  results  of  the  experi¬ 
mental  program  must  be  utilized  to  predict  the  required  wrap  angle  for  a 
desired  shock  front  velocity. 

Figure  I-XO  presents  a  curve  of  shock  front  velocity  versus 
wrap  angle,  this  curve  was  obtained  from  a  least  squared  fit  of  the  four 
data  points  which  were  available.  There  was  a  great  deal  of  scatter  is  the 
experimental  results  and  the  curve  is  only  accurate  to  *^Q0  ft/sec.  According 
to  figure  I- 10 ,  a  wrap  angle  of  9.2  degrees  will  produce  a  shock  front 
velocity  of  L,0OO  ft /sec. 
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e.  Overburden  Caleulalig&a 

A  pressure  pulse  of  extremely  abort  duration  would  fee  produced  if 
the  prlsaaeord  matrix  were  exploded  is  the  open  sir,  However,  the  total 
duration  can  be  sign!  ficaotiy  increased  and  the  major  portion  of  the  decay 
behind  the  frsast  ©oat rolled  if  an  overburden  is  utilised.  The  overburden 
will  serve  upward  after  the  detcc&tioo  takes  place  and  the  explosion  products 
will  expand  1st©  the  increasing  volume.  The  initial  volume  of  tins  eaplfe&ios 
cavity ,  the  peak  overpressure,  the  eversurden  mass ,  sad  the  expsr.sl  or,  proe»s» 
all  coat  rcl  the  aoveaeot  of  the  overburden.  A  mamaer  of  sistpi.ifyi&g,  assump¬ 
tions  ss*st  fee  made  before  the  sevMmnt  of  the  overburden  cam  fee  calculated . 

Figure  I -11  illustrates  the  e»4»i  which  was  used  for  the  overburden 
calculations .  A  time-varying  pressure  is  cos  fi  ted  is  a  cavity  of  Its  Halt# 
lateral  extent  which  has  ac  initial  height  of  *  .  The  lover  boundary  ©f 

V 

the  cavity  is  assumed  to  be  a  rigid  surface.  The  upper  boundary  of  the 
cavity,  which  is  formed  ay  the  overburden ,  is  ass  used  to  move  only  io  the 
vertical  direction  with  rigid  body  motive.  The  overburden  is  defined  fey  a 
©as a  per  unit  area,  a.  Only  one  coordinate  is  required  to  define  the  rigid 
body  motion  of  the  overburden.  Thus: 

x(t)  *  dlsplaeeaeat 
d2 

— *  acceleration 
dt 

The  expansioe  of  the  gases  in  the  cavity  can  be  described  fey 

r*  * 

P(t)  «  fp  ♦  al  1 - 

'  °  j  (Xo) 


where 


?o  *  peak  shock  overpressure 
a  *  atmospheric  pressure 
x  *  displacement 

»  initial  cavity  height 
y  »  ratio  of  specific  heats  *  1.3 
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Using.  Mewtoe’s  second  law. 


f  *  i  ^-r 
dt2 


the  following  equation  ssc  be  derived: 


A 


fp  ♦  *1 

1 

1  0 

|  m 

f 

*  *  «r.l 


♦ 


£  *  0 


this  nonlinear  differential  equation  was  solved  by  numrlcal  integration. 

figure  1-12  presents  the  results  of  the  overburden  calculations  for 
an  initial  overpressure  of  300  psi.  The  time  to  one-half  peak  pressure  has 
been  plotted  as  a  function  of  the  two  variables:  overburden  weight  tnd 
initial  cavity  height.  It  can  be  sees  that  increasing  either  the  initial 
cavity  height  or  the  overburden  weight  will  increase  the  tiae  to  one-half 
pressure. 

Results  from  Phase  1  indicate  that  the  duration  to  one-half  pressure 
is  appro xi oetely  fcl  percent  of  the  theoretical  value.  However,  the  effect  of 
expert  went  size  on  this  factor  was  not  known,  since  all  of  the  previews  tests 
had  been  conducted  in  approximately  the  saae  size  facility.  It  was  postu¬ 
lated  that  edge  effects  played  an  important  role  in  the  Phase  I  experiments, 
and  cm  the  basis  of  judgment,  a  correction  factor  of  90  percent  was  selected 
rat  her  than  the  experimental  value  of  percent.  Figure  1-12  was  then  used 

to  select  an  initial  cavity  height  of  36  inches  and  an  overburden  weight  of 

2 

500  lbs /ft  to  give  an  expected  time  to  one-half  pressure  of  37  msec. 

Figure  1-13  presents  the  theoretical  overpressure-t ine  history 
calculated  for  the  selected  parameters.  The  results  for  a  l-*ff  nuclear 
weapon  have  also  been  plotted  on  the  saae  figure  for  comparison ,  It  can  be 
seen  that  the  ^riaacord  curve  matches  the  nuclear  curve  only  in  the  neighbor¬ 
hood  of  the  one-half  peak  pressure.  However,  it  must  be  realised  that  the 
theoretical  priaacord  curve  accounts  only  for  the  pressure  decay  caused  by 
overburden  anti  on,  and  the  interaction  of  shocks  at  the  front  will  cause  a 
decay  immediately  behind  the  shock,  which  is  not  considered  in  these  eaiuc- 
latiocs,  dtier.  this  is  taken  into  account,  the  general  characteristics  of 
the  two  curves  should  be  quite  similar  down  to  one-half  peak  pressure.  The 
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pritta&ord  ourvo  teas  a  total  &ir«tios  of  160  ssao  acd  a  total  ispalse  of  15 


p% i  -sec. 
the  first 


fhs  auclear  cur»  has  a  total,  impulse  of  32  ps  i  -sec  j  hove*«r  *  in 
160  assse,  it  teas  an  impulse  of  appro xiaately  i@  ns  i  -sec . 
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Figure  I-;?.  Ai r* ftt ml- Loading  from  a  1  MU’  Surface  Burst 
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Figure  1-5.  Rarefaction  Waves 


Fiftur*  1-6.  Stres#  Wive#  in  «n  Eiwtlc  Hlif 


Figure  I-?*  Air-SIesi-Leduced  Groms  d  Mottoes  in  Superseisaic  Region 
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DEStREO  OmCCTlON  OF  TRAVEL 


Figure  1-8.  Pri**cord  Matrix 
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Figure  I *>9.  Peak  Shock  Overpr«#aur«  versus  Dce4  Density 
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Figure  I-li.  Hodei  for  Overburden  Calculations 
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Figure  1-12,  Sesults  of  Overfeuraets  Celcualtioas 
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Ki«ur*  1-13.  CoBpari »on  of  Pressuiw-Ti am  HiatoHee  for  Prliwujord  and  Nuclear 
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1.  Operations 

at.  Setae afcia®  Facility  Construction 


The  detessstion  facility  consists  of  si-eel  &feeet  piling  end  earth 
retaining  vails ,  a  wood  and.  steel  support  structure*  and  sand  overburden. . 
Figure  II -i  is  a  general  plan  view.  Figures  II -2,  II- 3,  and  11-1  are  three 
typical  sections  which  illustrate  the  staler  coaqx»ents  of  the  facility. 

Ihe  coast  ruction  of  the  detonation  facility  was  aeeoaplisaed  in  the 
fallowing  major  steps : 

Cl)  The  test  area  was  excavated  to  the  required  ddaeRsions. 

(2)  Steel  sheet  piling  vas  drives  around  the  perimeter  of  the 
facility,  except  for  a  small  access  roadway. 

(3)  Steel  H  piles  were  driven ,  and  tfF  beaas ,  to  span  the  model 
section, were  positioned  and  welded  in  place. 

ik)  lie  rows  of  prefab  ideated  footings,  ccluaass ,  and  girders 
were  set  into  place  and  nailed  together. 

(5)  The  remaining  steel  sheet  piling  was  driven  to  close  the 

roadway - 

{6}  The  pr inacord  raeits  were  installed,  tied  together,  and 
connected  to  the  plane  wave  detonation  generator. 


(I)  ’She  partially  prefabricated  wood  panels  were  positioned  an 
top  of  the  girders,  nailed  in  place,  and  filled  with  sand;  and  the  top  ply¬ 
wood  aeabers  were  stapled  to  the  floor  Joists. 

(%)  Sand  fill  was  placed  between  the  steel  sheet  piling  and  the 
wooden  support  structure. 

t'5)  Finally,  the  sad  overburden  was  placed  siasiitaneously  vi tn 
the  earth  bade  fill  which  surrounds  the  facility. 


The  facility  vas  constructed  in  a  rectangular  excavation  approximately 
120  feet  by  170  feet  in  plan  and  10  feet  deep.  As  originally  conceived,  the 
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facility  would  have  been  constructed  on  a  level  area  wills  the  earth  retaining 
well  is  &  bermed  ec®  figuration.  However,,  a  IS- foot  eat  was  required  for  this 
test  in  order  to  place  the  floor  of  the  facility  at  mb  elevation  below  rock 
layers  whies  were  present  at  the  site. 

The  steal  sheet  piling  was  drives  around  the  periaeier  of  the  test 
area  with  a  6-iheb  space  between  the  support  structure  arc  the  iaside  face 
of  the  sheet  piling.  TSsis  6-i»cfc  space  was  Isaac  filled  with  seed  to  provide 
additional  lateral  stability  to  the  facility.  She  steel  sheet  piling  was 
drives  to  a  depth  of  5  feet  and  was  cat  off  at  a  height  of  10  1/2  feet .  This 
was  to  provide  a  vertical  wort  log  face  and  to  support  the  unbalanced  hori  roots! 
force  which  results  frost  the  different  heights  between  the  overburden  and  the 
backfill. 

The  earth  retaining  wail  is  required  to  prevent  side  venting  after 
the  detonation  takes  place.  Side  venting  would  release  energy  fro*  the 
experiment ,  distort  the  oae -dimensional  awtioo  of  the  overburden,  and  present 
additional  safety  hazards.  The  earth  retaining  wall  was  constructed  by  back¬ 
filling  between  the  sheet  piling  and  the  side  of  the  opes  cut.  The  backfill 
material  was  placed  uncaapaetcd.-  It  consisted  of  a  predominately  silty 
sate  rial  which  was  present  at  the  site.  T&e  backfill  operation  was-  accoep- 
lished  sis&slta&ecusly  with  the  placement  of  the  sand  overburden  to  avoid 
excessive  bending  of  the  sheet  piling. 

The  support  structure  serves  three  purposes  simultaneously:  it 
provides  the  ini tied  cavity  in  which  the  explosion  takes  place.  It  positions 
the  priaecord  in  a  horizontal  plane  above  the  ground  surface ,  and  it  supports 
the  overburden.  The  initial  cavity  height,  free  the  surface  of  the  ground 
to  the  bottom  of  the  floor  panels,  is  3  feet.  After  the  detonation  takes 
place,  the  explosion  products  fill  this  initial  cavity  and  create  the  required 
peak  overpressure . 

Figures  I I -2  and  IX-3  are  typical  sections  which  illustrate  the 
construction  of  the  wood  port  on  of  the  support  structure.  The  bays  of 
footings,  eoiusrs ,  and  girders  were  placed  in  rows  which  span  the  length 
of  the  facility.  The  exterior  bays  were  located  on  k-foot  centers  while 
the  interior  bays  were  located  on  6-foot  centers.  The  columns  were  placed 
on  3- foot  centers,  along  the  length  of  each  row. 
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Sts  tvo  rows,  whs  ch  smks  up  the  central  bay,  were  broken  across  & 

'20- foot  span  at  the  Model  location,  and  steel  construction  was  used  is  place 
of  the  wood.  Steel  8  piles  were  drive©  to  depths  of  appro  xiaaiely  35  feet 
at  the  four  comers ,  a/ii  steel  *?  bstam  were  used  to  spaa  the  distance  be twees 
the  piles  is  earn  row.  Figure  tl-A  illustrates  this  portico  of  the  construc¬ 
tion.  Use  overburden  weight  is  transferred  to  the  M  piles  which  distribute 
the  load  along  their  length  instead  of  to  the  surface  of  the  great id.  this 
avoids  the  ©reloading  of  the  upper  portico  of  the  structural  aodel .  including 
the  door. 

the  wood  panels »  which  span  the  rows  of  girders*  were  partially 
prefabricated  prior  to  placement;  i.e. .  the  bottoe  plywood  sheets  were  stapled 
to  the  wood  Joists.  After  the  panels  were  sailed  into  pcs  it  ice  is  the 
facility,  the  space  between  the  Joists  was  filled  with  a  dense  cwaaerciai 
coa crete  sand  at  a  density  of  approximate ly  105  Is/ ft  .  The  dense  sand 
controls  the  initial  act  resect  of  the  overburden  by  prevent  lag  coapactioo  of 
the  lower  surface.  Shis  is  to  ensure  that  rigid  body  actions  ,  which  have 
been  assumed  fpr  the  overburden  calculations,  actually  take  place.  After 
the  sand  was  leveled,  the  top  plywood  sheets  were  stapled  in  place. 

flat  sand,  for  the  remainder  of  the  overburden  was  obtained  from  a 
local  borrow  pit.  It  was  selected  on  the  basis  of  obtaining  an  «n compacted 
density  of  approziaately  100  lb/ ft"  or  greater.  Use  sigh -density  require¬ 
ments  are  again  to  ensure  that  rigid  body  actions  take  place  and  to  produce 
the  ainimua  height  for  a  given  weight  of  overburden .  the  sane  overburden 
was  placed  from  the  exterior  of  the  test  area  utilising  a  long  booaed  crane. 

The  final  leveling  was  accomplished  by  a  saall  tractor  and  hand  labor.  The 
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total  weight  of  overburden ,  including  the  panel  sash,  was  500+5  lb/ ft  . 

Figures  II -5  through  Il-T  show  the  construction  sequence.  Figure  1 1 -5 
shews  the  detonation  facility  in  an  early  stage  of  the  construction.  The 
sheet  piling  arcus-s  the  perimeter  of  the  facility  can  be  seen  in  the  back¬ 
ground.  The  workmen  are  trenching  for  the  footings,  and  portions  of  the 
coluana  are  already  it  place.  Figure  II-c  shows  the  upper  portion  of  one  of 
the  a  piles  and  the  *ST  bean  which  spans  the  acdel  section.  Figure  II-T 
scows  che  completed  rows  of  footings,  columns  and  girders  .  The  primaeord 
racks  are  being  installed  before  placement  of  the  partially  prefabricated 


AWL  TO-65-26,  VOi  1 


wood  panels .  Use sc  panels  can  be  seen  stance 4  is  the  background, 
b.  Priaaoord  Installation 

i?5-gr*is ,  piast.i e—eoated  primaeord  manufactured  by  the  &asign~ 
Bickford  Ccapaey  ires  ati  iiz«4  for  this  tee t.  Teste  were  conducted  or.  300- foot 
eleaests  of  this  prisa.c©rd  and  it  wan  dgSarsaLaed  that  the*  datoe sties  velocity 
was  21  *000  ft /see  +1  percent. 

Cl)  Frimacord  Hack  Fabrication 

Ube  prisaccrd  matrix  parameters  were  selected  from,  the  data 
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presented  in  section  I.  A  lead  density  of  0.0?o  lb/ ft  and  &  wrap  «KSgle,  8, 
of  $.2  degrees  were  selected  to  respectively  produce  a  peak  equilibria*, 
overpressure  of  JOS'  psi  and  a  shock  front  velocity  of  k,®00  ft /sec. 


ft/sec. 


The  prissaecrd  sad  rices  were  fabricated  so  ■* seden  racks.  Use 
racks  were  constructed  of  2-incn  fey  %-inch  lumber  with  22-gage  metal  braces 
at  all  Joists  and  corners  (see  figure  11-3).  Tvc  rad  sires  were  required 
to  match  the  girder  spacing*  which  were  utilised  is  the  const  ruction  of  the 
detonation  facility.  IS*  larger  racks,  type  A,  were  ?  feet  2  inches  by 
21  feet  10  inches  and  the  smaller  racks ,  type  S,  were  2  fset  9  inches  by 
2**  feet  10  indies. 

Notches  were  cut  around  the  entire  perimeter  of  the  racks  to 
position  the  primacord.  Use  notch  spacing  was  determined  on  a  trial  and 
error  basis  in  an  attempt  to  simultaneously  match  the  two  parameters  of  load 
density  and  wrap  angle.  Table  II-l  summarises  the  notch  s pacings  which  were 
utilized.  Shis  notch  spacing  produced  a  wrap  angle  of  6.5  degrees. 


Rack  type 


Table  11-1 

SOTCH  SPACUtG 

Side  spacing 
{ inches ) 


a  spayc-ing 
Cinches) 


2.625 

2.5^ 
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The  wrapping  of  the  prijaeeord  was  accomplished  with  the  pria»eor4 
r*cks  s*g*ported  at  their  center  Use  with  *  chain  hoist,  Hsis  provided  a 
clear warn  for  free  passage  ©f  the  priasccrd.  The  priaacord  was  initially 
secured  to  *  top  coraer  of  the  rack  with  tape  «sd  stretched  to  the  opposite 
side  (bot&oa)  of  the  rack  m 4  placed  la  the  proper  acids,  ffee  cord  was  then 
passed  under  the  rack  and  again  stretched  to  the  opposite  side  (top)  and 
placed  in  the  proper  notch.  The  ecrd  was  alternately  passed  under  and 
thrown  over  the  top  of  the  rack  to  enable  the  wrappers  to  use  a  coat  mucus 
strand  of  priaacord.  Each  succeeding  wrap  cycle  was  a  repeat  of  the  previous 
operation  until  the  center  line  of  the  rack  was  reached.  At  this  point , 
the  person  throwing  the  cord  over  the  rack  proceeded  to  the  opposite  end  or 
the  rack  where  he  continued  Ms  operations .  This  procedure  continued  until 
the  entire  rack  had  been  wrapped  with  prisacord.  Figure  II-9  is  a  sketch  of 
this  operation. 


One  hundred  eighteen  thousand  feet  of  primacord  were  used  in 

fabricating  the  racks  and  lacing  the*  together  after  they  were  positioned 

in  the  facility .  The  total  volume  of  the  explosion  cavity  minus  the  voluae 
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of  the  wood  in  the  cavity  was  10,100  ft  .  The  actual  loading  density  was 

(116,000  n)(l?5  ®“~]  .. 

*  *  - — r-71 - *  0.0125  }~r 

j  7,000  ft'j  ft*' 

(2)  Priaacord  Rack  Placement 


The  racks  were  wrapped  with  priaacord  in  advance  of  their 
placement  in  the  detonation  facility.  Xkirisg  the  interis  period,  they  were 
stored  in  an  off-site  explosive  ordnance  bunker.  Transportation  to  the 
site  and  installation  in  the  detonation  facility  were  accomplished  in 
consecutive  steps.  Scan  of  the  prin&eord  racks  were  loaded  cn  a  flatbed 
trailer,  transported  tc  the  site,  and  unloaded  onto  pallets.  The  pallets 
were  lifted  into  the  detonation  facility  by  a  cran*  and  the  racks  were 
placed  by  workaen  while  the  flatbed  trailer  returned  for  another  load. 

Figure  1 1— T  shows  the  placement  o?  a  rack,  ihe  tops  of  the  columns  provided 
a  ledge  vhicn  positioned  the  prmacord  racks  in  a  nrri zoetal  plane  2  feet 
above  the  ground  surface. 
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Ttee  priaacard  os  the  end  of  each  rack  was  lace  a  to  Its  as&te  on 
the  es-d  sf  the  adjacent  rack  with  l?5-grais  prisno&rd.  Is  addition,  tms- 
verse  ties  vers  aaoe  with  short  lengths  of  lT5~grain  pris&cord.  These  ties 
were  aude  between  the  prisasord  <xa  the  sides  ©f  awsj&ees&t  reeks  at  3 -foot 
intervals .  Figure  11-40  shows  the  tying  operation. 

C  3 )  Pri  s&co  rd  Betsts  at i  ©a 

4  96-fbst.  strand  of  ^-frai®  pries&cord  was  osaseeted  to  <@asfc 
individual  piece  of  prisacord  across  the  entire  width  of  the  mtrtx  ©a  the 
detonation  end  of  the  facility.  Forty-eight  strands  of  175 -krais  primteoni 
were  connected  to  the  lOO-gcais  prisacord  os  2-ftoot  centers.  Each  of  the 
«8  strands  was  placed  is  a  plastic  casing  which  passed  through  the  elywood 
bulkhead,  sheet  piling,  and  backfill  sate  rial  to  the  original  ground  surface 
(see  figure  IX-ilk 

Figure  11-12  illustrates  the  geossetric  patters  of  e^uai-lesutth 
115-graic  prissacord  strands  which  were  utilized  to  fabricate  the  plane -wave 
generator.  Delicate  blasting  caps  were  used  to  initiate  detonation.  Tfee 
detonation  wave  traveled  down  the  plane  wave  generator  and  simultaneously 
detonated  the  %00-grain  priaaccrd  at  hB  separate  points.  Phase  1  results 
Indicate  that  this  detonation  ache*  will  generate  *  wave  which  is  plane 
across  the  width  of  the  facility  within  a  few  si crose coeds . 

c.  Betosation  Sequence 

A  dry -rur,  of  the  countdown  procedure  *«#*£  conducted  before  the  be*?t 
to  ensure  that  all  instrument at iec ,  hissing,  and  firing  circuits  were  properly 
functioning.  The  countdown  procedure  consisted  of  the  following  steps: 

Cl)  At  T-l  test  area  was  cleared  of  all  personnel. 

(?)  The  blasting  caps  were  installed  on  the  plane  wave  generator 
and  tae  firing  lines  were  connected  to  Use  blasting  caps. 

(3)  The  firing  lines  were  connected  to  the  safety  box  at  the  c cowans 
center.  Before  the  firing  circuit  could  be  energised  fro*  the  power  supply, 
a  plug  ssui  to  be  paysicaliy  inserted  into  the  safety  box  and  a  button 
depressed  and  held  down. 


(*)  At  T-15  aiautes,  the  area  was  visually  checked  by  guards  to 
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ensure  that  the  are*  was  still  clear  of  personae* . 

C5)  At  f-5  minutes ,  all  Instrumentation  circuits  were  energised. 

(b)  At  T-2  minutes ,  the  firing  system  programmer  vas  initiated, 

C?3  At  f-l  w&mit e»  &  siren  van  sounded.  ©se  ping  was  inserted  late 
tide  safety  box  and  the  button  was  depressed. 

CS1  At  T-I5  seconds ,  a.  verbal,  countdown  vas  started. 

(9)  At  ?~1Q  seconds,  the  assual  calibration  switch  was  thrown. 

CIO}  At  ?-8  seconds,  the  call  brat  lor.  switch  vas  opened. 

Ul)  At  f~-2  seconds,  the  cameras  were  turned  on. 

(12}  At  f-0  seconds ,  the  explosives  were  detonated. 

(13)  'She  facility  vas  visually  inspected  for  us  detonated  explosives 
and  the  “All  Clear”  siren  vas  sounded. 

d.  Safety 

The  command  center  vas  located  approximately  >00  feet  from  the  edge 
of  the  detonation  facility.  Only  essential  personnel  were  allowed  in  the 
command  center.  Spectators  were  not  allowed  to  approach  within  1,000  feet 
of  the  betas siioc  facility  frost  T-I  hours  until  after  the  “All  Clear5*  siren 
was  sounded, 

The  weather  criteria  were  based  on  a  maxi  aum  vied  velocity  of  10 
siles  per  hour  with  a  heading  such  that  sand  or  debris  would  set  be  blown 
toward  an  active  portion  of  the  Sirtland  Air  Force  Base. 

Because  of  the  siting  on  an  Air  Force  base,  close  coordination  was 
required  with  the  FAA  tower  to  ensure  that  aircraft  did  not  pass  over  the 
facility  in  their  landing  and  take  off  patterns  during  the  firing  sequence. 

2.  Instrumentation 

a.  General  Layout 

Figure  11-13  illustrates  the  general  instrumentation  layout.  The 
locations  of  the  structural  model  and  the  hardened  cable  project  are  indicated. 
The  general  location  of  the  air  pressure  gages  .  the  sear-surface  soil  pressure 
gages,  and  the  drill  holes  for  the  earth  free- field  gages  are  also  shown. 

i<6 
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The  specific  instrumentation  details  are  reported  Is  Paris  2  and  3.  A 
detailed  instnsaeiitaiion  list  is  provided  is  appendix  A. 

All  the  instrumentation  cosies  were  brought  oat  of  the  detoB&tie® 
facility  ia  five  aajsr  t reaches  wfeiefe  has  a  sd&isgasB  depth  of  cover  ©f  3  feet 
of  earth.  The  sheet  piling  immediately  over  the  cable  t rescues  woo  drives 
to  1— foot  depth;  is® teaus  of  the  ssem*!  5  feet  ts  prerest  shearing  of  the 
cables  at  this  biaeosiisuity  -  The  cables  were  farther  protected  fey  he  lag 
placed  is  plastic  ena&ults  wfeiefe  ran  ©OKtims^aaly  frost  all  laatraiaeat-afei as 
locations  to  the  exterior  of  the  deioeatian  facility,  the  cables  were  brought 
to  the  surface  approximately  $0  feet  frtas  the  north  edge  of  the  facility  and 
froc  this  point  they  ran  above  ground  to  the  iosiroaest ahi on  vans . 

b.  Instrumentation  Recording 

■23®  electronic  signals  produced  by  the  free-fieM  gags  a  is  this 
experiment  were  recorded  is  the  Air  Force  Weapons  laboratory  Field  lastm- 
meniatioa  Trailer  which  was  located  app  roximate ly  5-00  feet  to  the  north  ©f 
the  test  site.  The  C3EC  “System  B*  <J~kc  carrier}  was  used  for  signal  csssdi- 
tloaisg  of  the  velocity  gages.  “She  remaining  free- field  gages  were  operated 
by  Allegheny  Instrument  "Sensor  Analogue  Modules'’'  {SAM-1}  signal  condi  t  lotting 
equipment.  The  SAM-1  is  a  DC  system  incorporating;  a  power  supply  and  a  SC 
wide  band  amplifier  for  each  gage.  Four  a apex  CP-100  and  two  {StC  Pfe-330© 
tape  recorders  were  used  to  record  all  free- field  instrumentation  with  the 
exception  of  tbe  earth  tiae  of  arrival  system ,  watch  was  recorded  fey  fast  ax 
camera.  A  time  code  generator  was  used  to  place  an  IBIS  time  code  ca  each 
tape.  Simulated  calibration  signals  were  sect  through  each  recording 
channel  lamed  ate  ly  before  and  after  the  shot  by  shunting  selected  resistors 
across  one  ara  of  the  transducer  bridge .  Figure  II-li  is  a  schematic  of  the 
ins  truaent  ati  oe  reccrditg  systea. 

c .  Camera  Coverage 

high-speed  notice  picture  ar.d  Pastax  cameras  were  utilised  to 
photograph  the  event,  these  cameras  were  mounted  or.  two  wooden  starts  ,  one 
to  the  north  and  one  to  the  west  of  the  detosaticn  facility,  Wooden  targets 
were  positioned  to  give  the  relative  movement  between  the  tap  and  Pottos  of 
the  overburden  as  illustrated  is  figure  11-15 .  the  camera  ©overage  wan 
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planned  to  give  data  on  the  overburden  motion  and  any  venting  of  the  explosion 
products  as  veil  as  to  provide  toted  pictorial  coverage  of  the  event. 
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Figure  II-l.  Plan  View  of  Detonation  Facility 
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Figure  II-2.  Cross  Section  of  Detonation  Facility 
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Figure  II-3.  Cross  Section  of  Detonation  Facility 
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Figure  II-4.  Cross  Section  of  Steel  Support  Structure 


Figure  II-5.  Early  Stage  of  the  Facility  Construction 
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Figure  II-J.  Primacord  Rack  Placement 
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Figure  II-9.  Primacord  Wrapping 
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SECTION  III 
RESULTS 

This  experiment  produced  a  peak  overpressure  of  312  psi ,  a  time  to  one- 
half  peak  pressure  of  18.2  msec  and  a  total  duration  of  170  msec.  A  composite 
pressure-time  history  curve  is  given  in  figure  IX-1  which  defines  the  air 
free  field  environment.  The  impulse  contained  under  this  curve  is  approximately 
10.1*  psi-sec.  The  average  shock  front  velocity  was  5*120  ft/sec  +11  percent, 
and  the  wave  front  was  essentially  plane,  across  the  central  portion,  as  it 
traveled  the  length  of  the  facility. 

The  overburden  moved  upward  with  rigid  body  motions  until  it  cleared 
the  top  of  the  metal  sheet  piling.  The  overburden  continued  to  move  upward 
as  a  mass  until  it  reached  a  height  of  approximately  125  feet  at  the  firing 
end  {see  figure  III-l).  At  this  point,  the  individual  particles  03f-sand 
began  to  break  up  and  a  dust  cloud  continued  to  rise  to  nearly  twice  this 
height.  The  overburden  displacement  versus  time  measurements  taken  from  the 
Fastax  camera  coverage  of  the  target  stands  are  presented  in  figure  III -2. 

No  differential  motions  between  the  top  and  the  bottom  of  the  overburden 
could  be  detected  in  the  target  stand  data. 

A  small  amount  of  sidewall  venting  was  first  detected  at  the  southeast 
comer  of  the  facility.  This  venting  progressed  across  the  detonation  end  of 
the  facility  and  then  proceeded  down  the  sidewalls  and  across  the  reflection 
end  of  the  facility.  Major  venting  did  not  occur  until  the  bottom  of  the 
overburden  had  cleared  the  top  of  the  sheet  piling. 
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Figure  III-2.  Overburden  Motion 
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SECTION  IV 
DISCISSION 

The  peak  overpressure  and  shock  front  velocity  curves  developed  during 
Phase  I  (figures  1-9  and  1-10)  have  been  corrected  to  include  the  Phase  II 
data.  These  corrections  were  quite  small.  Figures  IV-1  and  IV-2  present  the 
corrected  peak  overpressure  versus  load  density,  and  shock  front  velocity 
versus  wrap  angle  prediction  curves. 

The  composite  air-pressure  curve  obtained  from  this  experiment  is  compared 
with  the  results  from  a  1-MT  nuclear  weapon  in  figure  IV-3.  It  can  be  seen 
that  the  desired  durations  were  not  obtained  in  this  experiment.  The  observed 
duration  to  one-half  peak  pressure  was  49.2  percent  of  the  desired  duration 
and  44.4  percent  of  the  computed  duration.  This  agrees  well  with  the  43.9 
percent  observed  in  the  Phase  I  experiments.  It  had  been  postulated  that 
edge  effects  played  a  major  role  in  producing  the  observed  low  efficiency 
factor.  However,  this  does  not  now  appear  to  be  the  ease  since  the  perimeter 
length  to  plan  area  ratio  was  decreased  from  0.15  for  the  Phase  I  facilities 
to  0.0342  for  the  Phase  II  facility  with  no  substantial  increase  in  the 
efficiency  factor. 

Figure  IV-4  presents  the  displacement-time  history  of  the  overburden 

s 

motion.  A  comparison  of  the  theoretical  and  observed  time  histories  reveals 
that  the  overburden  was  initially  accelerated  much  more  rapidly  than  was 
computed,  while  at  late  times  there  is  fair  agreement.  The  observed  high 
initial  acceleration  will  cause  a  much  faster  decay  of  pressure  than  was 
computed. 

Since  the  duraticaa  efficiency  factors  from  both  Phase  I  and  Phase  II 
agree  and  experimental  results  clearly  indicate  that  the  overburden  calcu¬ 
lations  are  in  error,  it  is  obvious  that  these  calculations  must  be  adjusted 
in  order  to  successfully  make  predictions  for  future  experiments.  Figure 
IV-5  has  been  extracted  from  reference  1.  It  is  a  plot  of  duration  efficiency 
factors  to  one-half  peak  pressure  for  various  cavity  heights.  The  results 
in  figure  1-12  have  been  corrected  to  include  the  efficiency  factors  given 
in  figure  IV-5.  The  corrected  set  of  curves  are  presented  in  figure  IV-6. 
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Figure  IV-1.  Peak  Shock  Overpressure  versus  Load  Density  (Corrected) 


PRIMACORD  WRAP  ANGLE, decree* 

Figure  IV-2.  Shock  Front  Velocity  versus  Wrap  Angle  (Corrected) 
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Figure  IV-3.  Comparison  of  Pressure-Time  histories  for  Primacora  and  Nuclear 
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SECTION  V 

CONCLUSIONS  AND  RECOMMENDATIONS 

The  parameters  of  the  simulation  technique  scaled  exactly  between  the 
small  Phase  I  tests  and  the  large  Phase  II  test.  The  peak  overpressure  and 
shock  front  velocity  prediction  curves,  which  were  developed  during  Phase  I, 
appear  to  be  valid  within  approximately  10  percent  and  only  small  corrections 
were  required  to  include  the  Phase  II  data. 

The  observed  durations  to  one-half  peak  pressure  were  only  l»9.2  percent 
of  those  desired.  Duration  efficiency  factors  have  been  applied  to  the 
theoretical  calculations  to  account  for  this  discrepancy.  Figure  V-l 
summtu'izes  these  calculations  and  illustrates  the  overburden  versus  initial 
cavity  height  required  to  simulate  a  300-psi  shock  wave  for  a  1-MT  detonation. 

For  future  300-psi  simulations ,  it  is  recommended  that  a  load  density  of 
0.0720  pound  of  PETN/ft  and  a  wrap  angle  of  8.6  degrees  be  used.  The  initial 
cavity  height  and  the  overburden  mass  should  be  selected  from  the  values 
given  in  figure  V-l.  The  primacord  matrix  should  be  positioned  2k  inches 
above  the  ground  surface  when  a  36-inch  cavity  height  is  utilized.  Similar 
reconmendations  cannot  be  made  for  larger  cavity  heights  until  additional 
tests  are  conducted. 


Figure  V-l.  Overburden  versus  Cavity  Depth  Required  to  Simulate 


This  page  intentionally  left  blank. 
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SECTION  VI 
INTRODUCTION 


1.  Objectives 

The  air  free-field  measurements  were  made  to  determine  the  air-blast 
pulse  produced  by  the  experiment,  to  provide  input  data  for  the  structural 
model  and  the  earth  free  field,  and  to  provide  data  to  improve  prediction 
techniques  for  future  tests.  In  addition,  measurements  were  made  of  the 
shock  wave  which  was  released  to  the  atmosphere. 

2.  Background 

The  research  performed  during  the  first  phase  of  this  project  was  used 
as  a  basis  to  design  this  large-scale  experiment.  During  Phase  I,  four 
primacord  tests  were  conducted  and  experimental  data  were  collected  to 
determine  the  load  density,  primacord  wrap  angle,  cavity  height,  and  mass 
of  overburden  which  are  required  to  produce  a  desired  simulation.  As  outlined 
in  section  I,  an  experimental  approach  was  necessary  since  the  analytical  . 
problem  is  nearly  intractable. 

Many  of  the  Phase  I  conclusions  were  arrived  at  with  a  minimum  of  data 
points.  Also,  all  of  the  previous  experiments  were  conducted  on  a  much 
smaller  scale  than  the  Phase  II  experiment.  This  experiment  will,  therefore, 
provide  valuable  data  for  improving  the  prediction  techniques  for  future 
simulations.  For  example,  it  was  found  in  Phase  I  that  a  36-inch  cavity 
would  produce  a  duration  which  was  UI4  percent  of  the  theoretically  computed 
duration.  The  effect  of  scaling  the  loaded  area  on  this  parameter  can  be 
determined  on  the  basis  of  the  Phase  II  results. 

3.  Predictions 

As  discussed  in  section  VI-2,  the  Phase  I  results  were  utilized  to 
predict  the  parameters  required  to  give  the  desired  simulation.  In  addition, 
judgment  was  utilized  to  determine  the  effect  of  increasing  the  scale  of 
the  experiment.  On  this  basis,  the  following  predictions  were  made: 

Peak  overpressure  300  psi 

Time  to  one-half  pressure  37  msec 
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Total  duration 
Shock  front  velocity 


160  msec 
l*, 800  ft/sec 


SECTION  VII 


PROCEDURE 


1.  Instrumentation 

a.  Air-Blast  Gages 

The  prims  instrumentation  for  this  portion  of  the  experiment  consist*’-! 
of  flush-mounted  pressure  transducers.  These  gages  provide  pressure-time 
histories  as  well  as  shock. front  time-of-arrival  data.  Two  types  of  pressure 
transducers  were  used,  the  Norwood  Model  111  and  the  Scheavitz-Bytrex  Model 
HFG-2000 .  The  Norwood  transducer  is  a  conventional  strain  gage  device 
selected  because  it  proved  to  be  dependable  and  rugged  during  the  Phase  I 
test  series.  The  main  weakness  of  this  gage  is  its  relatively  low  frequency 
response  (20,000  cps)  which  causes  a  moderate  amount  of  ringing  when  the  gage 
is  subjected  to  an  air  shock.  The  Scheavitz-Bytrex  gage  is  a  recently 
developed  pressure  transducer  which  uses  a  semiconductor  strain  gage  to  link 
a  very  stiff  diaphragm  to  the  transducer  body.  This  gage  was  selected 
because  of  its  relatively  high  frequency  response  (50,000  cps).  Calibration 
was  accomplished  in  the  field  with  static  air  pressure. 

b.  Shock  Front  Time  of  Arrival 

A  line  of  break  wires  was  placed  along  the  length  of  the  facility  to 
provide  additional  shock  front  time-of-arrival  data.  The  break  wires  were 
manufactured  in  the  Air  Force  Weapons  Laboratory  from  very  lightweight  foil. 
The  break  wires  wer^  connected  to  a  DC  circuit  such  that  when  a  foil  element 
was  broken,  a  step  increase  in  voltage  was  produced.  The  seven  break  wires 
were  multiplexed  into  a  single  channel  for  recordings 

c .  Mi crob  aragraphs 

Microbaragraphs  were  used  to  measure  the  intensity  of  the  shock  wave 
formed  in  the  surrounding  air  when  the  facility  vents.  These  self-contained 
instruments  had  a  variable-- reluctance  bridge  pressure  sensor  head  supplied 
by  a  1,000-cps  carrier,  preamplifier,  detector,  DC  amplifiers,  and  Brash 
recorders.  The  sensing  head  is  a  Wianko  twisted  Bourdon  tube. 
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2,  Gage  Placement 


The  pressure  transducers  were  mounted  in  concrete  blocks  which  had  a 
concave  upper  surface  to  facilitate  the  attachment  of  calibration  fittings. 
-The  blocks  were  set  in  the  field  with  the  top  of  the  gage  flush  with  the 
ground  surface.  The  gage  leads  were  threaded  through  a  flexible  metal,  tubing 
which  had  been  cast  in  the  concrete.  After  exit  from  the  concrete,  the  cables 
were  protected  by  plastic  conduit.  The  upper  portion  of  the  mounting  block 
was  grouted  level  with  the  ground  surface  when  the  calibration  of  each  gage 
was  completed.  See  figure  VII-1. 


Locations 


Figure  VII-2  is  a  plan  view  of  the  test  area  showing  the  location  of  the 
air  free-field  instrumentation.  It  will  be  noted  that  each  position  containin 
a  pressure  transducer  has  been  assigned  a  location  number.  The  actual 
distance  between  various  gage  locations  was  surveyed,  and  these  distances 
have  been  indicated  in  figure  VII-2.  The  break  wires  were  placed  at  25-foot 
intervals,  commencing  at  the  detonation  end  of  the  facility.  Their  location 
will  be  designated  by  the  letters  D.  E. ,  plus1  a  distance  to  the  gage  of  intere 
i.e,,  D.E.  +  75  feet.  Table  VII-1  lists  the  transducers  used  at  each  of  the 
locations  shown  in  figure  VII-2. 


The  microbaragraphs  were  located  along  a  line  at  distances  of  UOO,  800, 
and  1,200  feet  from  the  north  edge  of  the  test  facility. 
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Table  VII-1 
INSTRUMENTATION  LIST 


Location  number 

Instrument  type 

1 

Norwood  Model  111 

1* 

Norwood  Model  111 

5 

Norwood  Model  111 

7 

Norwood  Model  111 

10 

Scheavitz-Bytrex  Model  HFG-2000 

14 

Norwood  Model  111 

18 

Norwood  Model  111 

22 

Scheavitz-Bytrex  Model  HFG-2000 

24 

Norwood  Model  111 

28 

Scheavitz-Bytrex  Model  HFG-2000 

30 

Norwood  Model  111 

32 

Scheavitz-Bytrex  Model  HFG-2000 

34 

Norwood  Model  111 

D.E. 

Break  Wire 

d.e.  +  25  ft 

Break  Wire 

d.e.  ♦  50  ft 

Break  Wire 

d.e.  +  75  ft 

Break  Wire 

d.e.  +  100  ft 

Break  Wire 

D.E.  +  125  ft 

Break  Wire 

D.E.  +  150  ft 

Break  Wire 
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Figure  VII-1.  Pressure  Transducer  Mount 
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SECTION  VIII 


RESULTS 


1.  Air  Pressure  Gages 

Nine  pressure-time  history  traces  were  recorded  from  the  13  air  pressure  I 

s 

gages  which  were  installed  in  the  experiment.  In  addition,  three  of  the  ; 

channels  {18,  22,  and  28)  contained  time-of-arrival  data,  but  pressure-time  f 

data  could  not  be  recognized  or  extracted  from  th_  excessive  amount  of  noise  ■ 

i 

which  was  recorded.  The  cause  of  this  excessive  noise  has  not  been  determined.  j 
The  gage  at  position  No.  30  was  disconnected  before  the  test  because  of 
faulty  signal-conditioning  equipment. 

Appendix  B  contains  the  time  history  records  for  all  of  the  air  pressure 
gages  which  recorded  usable  data.  Table  VIII-1  is  a  summary  of  the  shock 

front  time  of  arrival  at  each  of  the  gage  locations.  i 

! 

Table  VIII-1  ! 

j 

SHOCK  FRONT  TIME  OF  ARRIVAL  AT  PRESSURE  TRANSDUCER  LOCATIONS  ! 


Location 

No. 


Time  of  arrival 
(msec) _ 


15.8 
20.  h 
2U.9 


15.8 

20.h 

25.3 
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2.  Break  Wires 

All  seven  of  the  break  wires  performed  satisfactorily  and  the  times  of 
arrival  recorded  at  each  location  are  summarized  in  table  VIII-2. 


Table  VIII-2 

SHOCK  FRONT  TIME  OF  ARRIVAL  AT  BREAK  WIRE  LOCATIONS 


Location 

Dist  from  Detonation  End  (ft) 

Detonation  End 
Detonation  End  +25 
Detonation  End  +  50 
Detonation  End  +  75 
Detonation  End  +  100 
Detonation  End  ♦  125 
Detonation  End  +  150 

3.  Microbar agraphs 


Time  of  Arrival 
(msec) 

0.0 
5. It 
10.8 

16.1 

21.2 
25.9 
30.3 


The  results  obtained  from  the  mi crobaragraphs  are  summarized  in  table 
VIII-3. 


Table  VIII-3 
MICROBARAGRAPH  RESULTS 


Range 
( feet ) 

Peak  positive  phase 
overpressure 
(microbars)  (psi) 

Peak  negative  phase 
overpress  \ire 
(microbars)  (psi) 

Positive 

duration 

(sec) 

Total 

duration 

(sec) 

L00 

8,U00 

0.123 

It,  800 

0.0705 

0.20 

0.53 

800 

3,900 

0.057!+ 

2,1+00 

0.035 

0.20 

0.53 

1,200 

3,360 

0.0495 

1,680 

0.025 

0.20 

0.53 
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SECTION  IX 
DISCUSSION 


1.  Pressure-Time  History 

A  review  of  the  pres sure -time  history  records  presented  in  appendix  B 
will  reveal  the  presence  of  many  small  shocks  during  the  early  portion  of 
the  record  (approximately  the  first  10  msec).  Originally,  it  was  planned  to 
ignore  these  minor  shocks  and  to  wait  for  thermal  equilibrium  before  defining 
the  peak  overpressure.  However,  the  shock  portion  of  the  curve  represents 
too  much  of  the  inqpulse  to  be  ignored.  Since  the  soil  tends  to  integrate 
or  smooth  out  the  pressure  pulse,  it  was  decided  to  substitute  a  smoothed 
curve  for  each  observed  pressure  pulse,  taking  care  to  preserve  impulse.  To 
arrive  at  this  smoothed  curve ,  the  early  portion  of  each  individual  record 
was  divided  into  a  number  of  equal  time  increments.  A  mean  pressure  value 
was  selected  for  each  time  increment  such  that  the  impulse  was  preserved. 
Finally,  a  smooth  curve  was  constructed  utilizing  the  mean  values,  and  it 
was  assumed  that  this  fit  represented  the  pressure  pulse  seen  at  a  particular 
gage. 

Figure  IX-1  is  a  composite  pressure-time  history  curve  derived  from  the 
individual  smoothed  curves  described  above.  The  peak  pressure  of  the  composite 
curve  was  determined  by  averaging  all  of  the  peak  pressures  from  the  smoothed 
curves.  To  determine  the  wave  shape  of  the  composite  curve,  the  smoothed 
curve  representing  each  pressure-time  history  was  normalized  to  the  average 
peak  pressure.  These  normalized  curves  were  then  averaged  together  to 
complete  the  composite  curve. 

Three  of  the  pressure  records  indicated  negative  pressures  at  relatively 
early  times.  These  curves  were  adjusted  by  linear  base  line  shifts  such 
that  they  reached  a  zero  pressure  at  a  time  which  agreed  with  the  zero 
pressure  from  the  composite  curve. 

Some  Justification  exists  for  the  rather  complicated  derivation  of  a 
composite  air  pressure  curve.  Figures  XIV- 2  and  XIV-3  compare  the  pressure- 
time  histories  seen  at  two  shallow  buried  soil  stress  gages  with  the  composite 
air  pressure  curve.  The  shocks  have  already  decreased  in  magnitude  and  the 
measured  stress  curve  appears  to  be  approaching  the  shape  of  the  composite  curve 
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1  2.  Shock  Front  Velocity 

f 

Shock  front  velocities  have  been  computed  from  the  shock  front  time-of- 
|  arrival  data  and  the  results  have  been  plotted  in  figures  IX-2  through  IX-5 » 

|  respectively,  for  the  north  gage  line,  central  gage  line,  south  gage  line 
|  and  break  wires.  These  plots  contain  interval  velocities  as  well  as  average 
|  velocities.  The  distance  between  gages  was  measured  to  the  nearest  0.01  foot; 

I  however,  the  physical  location  of  the  gages  in  the  test  facility  was  only 

i 

|  known  to  +1  foot.  A  summary  of  the  interval  shock  front  velocities  is 
'  presented  in  figure  IX-6  along  with  an  average  shock  front  velocity  of  5,120 
ft/sec.  It  should  be  noted  that  there  appears  to  be  an  increase  in  velocity 
with  the  distance  traveled.  The  individual  values  agree  with  the  average 
value  within  +l6  percent.  However,  the  central  line  of  gages  represents'  the 
largest  deviation  from  the  average  value  in  both  the  plus  and  minus  direction. 
These  points  are  suspected  to  reflect  surveying  errors  since  their  average 

■  velocity  is  exactly  5,120  ft/sec.  If  these  points  are  disregarded,  the  indi¬ 
vidual  values  agree  with  the  average  value  within  +?1  percent. 

Figure  IX-7  is  a  plot  of  equal  time  contours  which  define  the  position 
of  the  shock  front  in  the  test  facility.  The  shock  front  is  essentially 
plane  across  the  central  50  feet  of  the  test  facility.  The  maximum  distance 

■  deviation,  for  a  given  time,  in  the  central  area  is  3  feet.  This  is  equiva- 
;  lent  to  approximately  0.6  msec.  There  is  a  definite  edge  effect  along  the 

;  entire  length  of  the  facility.  This  may  be  due  to  the  smaller  Drimaeord  racks 

! 

I  which  were  utilized  at  the  edge  positions. 

| 

1  3.  End  Wall  Reflections 

1 

|  When  the  shock  front  reaches  the  end  of  the  facility  a  reflected  shock  is 
|  formed  which  travels  in  a  direction  opposite  to  the  initial  flow  and  causes 
|  an  increase  in  pressure.  Table  IX-1  contains  a  summary  of  the  observed  results. 
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Table  IX-1 
REFLECTION  RESULTS 


where  t  *  time  of  reflected  shock  arrival 
r 

=  strength  of  reflected  shock 

The  reflected  shock  pressures  have  been  plotted  as  a  function  of  travel 
distance  in  figure  IX-8.  If  a  linear  fit  to  the  data  is  assumed,  a  reflection 
factor  of  approximately  2  is  obtained  at  the  end  wall  (reflection  factor  * 
incident  pressure  plus  reflected  pressure  divided  by  incident  pressure).  For 
a  shock  wave,  the  reflection  factor  from  a  perfectly  rigid  wall  is  6.5. 
However,  the  detonation  wave  is  not  a  true  shock  and  the  end  wall  provides  an 
extremely  soft  boundary.  Figure  IX-8  also  indicates  that  the  reflected 
pressure  should  have  died  out  by  the  time  it  has  traveled  60  feet.  No 
reflected  pressures  were  observed  on  the  pressure  records  at  70  feet  from 
the  reflection  end.  The  average  reflected  shock  front  velocity  has  been 
obtained  from  time-of-arrival  data  in  figure  IX-9.  The  average  velocity  of  . 

1* ,000  ft/sec  corresponds  to  a  pressure  of  200  psi  in  air,  which  agrees  fairly 
well  with  average  measured  reflected  pressures. 

It.  Shock  Pressure  Vented  to  the  Atmosphere 

In  figure  IX-10  the  air  shock  which  was  created  by  the  venting  into  the 
atmosphere  is  compared  with  an  equivalent  unconfined  hemispherical  high 
explosive  detonation  with  a  total  charge  weight  of  3,000  pounds.  It  can  be 
seen  that  the  Long-Duration  High -Explosive  Simulation  Technique  creates  an 
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I  air  shock  which  is  at  least  an  order  of  magnitude  lower  in  peak  overpressure 
1  than  its  equivalent  charge  weight  of  high  explosives . 
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Figure  IX-2.  Shock  Front  Velocities — Worth  Gage  Line 
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Figure  IX-3.  Shock  Front  Velocities — Central  Gage  Line 
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Figure  IX-6.  Summary  of  Shock  Front  Velocities 


Figure  IX-7.  Shock  Front  Positions 
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Figure  IX-8.  Attenuation  of  Reflected  Pressure 
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Figure  IX-9.  Reflected  Shock  Front  Velocity 
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SECTION  X 

CONCLUSIONS  AND  RECOMMENDATIONS 


The  composite  pres sure -time  history  curve  given  in  figure  IX-1  defines 
the  air  free-field  environment  which  was  produced  by  this  experiment.  It 
indicates  a  peak  overpressure  of  312  psi,  a  time  to  one-half  peak  pressure  of 
18.2  msec,  and  a  total  duration  of  170  msec.  The  impulse  contained  under 
this  curve  is  approximately  10. 1*  psi-sec. 

The  shock  front  velocity  produced  in  this  experiment  increased  with  travel 
distance  and  ranged  from  U,550  ft/sec  to  5,680  ft/sec.  The  average  velocity 
was  5,120  ft/sec  +11  percent.  The  wave  front  was  essentially  plane  in  the 
central  50-foot  portion  of  the  loaded  area. 

The  reflected  pressure  at  the  end  wall  is  of  the  same  order  or  larger 
than  the  incident  pressure.  The  reflected  pressure  attenuates  rapidly  with 
distance,  and  for  this  experiment  it  had  decayed  to  zero  in  60  to  70  feet  of 
travel. 

The  air  shock  which  was  created  in  the  atmosphere  when  the  experiment 
vented  was  an  order  of  magnitude  lower  than  that  created  by  an  equivalent 
weight  of  high  explosives  detonated  unconfined  in  a  hemispherical  charge. 


V/ 
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PART  3 


EARTH  FREE  FIELD 
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SECTION  XI 


INTRODUCTION 


1.  Objectives 

The  earth  free-field  measurements  were  made  to  determine  the  response 
of  a  known  medium  to  a  controlled  input,  to  determine  the  attenuation  of 
free-field  quantities  with  depth,  to  provide  input  data  for  the  structural 
model,  and  to  provide  data  for  determining  the  mass  of  material  realistically 
loaded  so  that  larger  follow-on  experiments  may  be  appropriately  scaled. 

2.  Background 

Earth  free-field  measurements  have  been  made  on  numerous  nuclear  and  high- 
explosive  detonations.  These  data  have  been  utilized  to  derive  a  number  of 
empirical  prediction  techniques  as  veil,  as  to  provide  check  points  to  verify 
analytical  solutions.  Any  theory  used  to  predict  the  motions  and  stresses 
induced  in  the  earth  from  an  air-blast  loading  requires  a  knowledge  of  the 
properties  of  the  soil.  Therefore,  a  comprehensive  site  exploration  program 
was  accomplished  to  determine  the  properties  of  the  soil  at  the  test  site. 

3.  Soil  Survey 

The  site  exploration  program  consisted  of  seismic  refraction  and  uphole 
velocity  surveys,  the  obtaining  of  undisturbed  samples  and  the  laboratory 
testing  of  samples.  Laboratory  investigations  included  mechanical  analysis , 
constrained  modulus  tests,  triaxial  tests,  and  moisture  content  and  density 
determinations . 

The  results  of  the  seismic  survey  are  summarized  in  table  XI -1  (reference 
5).  The  nature  of  higher  velocity  interface  detected  at  268  feet  is  unknown 
but  it  is  probably  the  surface  of  an  older,  denser  alluvial  material.  No 
higher  velocities  were  detected  to  a  depth  of  1,500  to  2,000  feet. 
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Table  XI-1 
SEISMIC  PROFILE 


0-58 

1700 

58-268 

2550 

268-2 ,000 

67^0 

The  teat  site  consists  of  layered  media  ranging  from  a  dense  sand  to 
medium  gravel  with  some  mixing  of  clay  and  caliche.  There  was  no  ground 
water  table  above  90  feet.  The  soil  profile  is  shown  'n  figure  XI -1  (reference 
6).  The  boring  denoted  by  "HE"  was  sampled  by  using  5-inch  Shelby  tubes  to 
obtain  1-foot  undisturbed  samples  at  5-foot  intervals  to  a  depth  of  approximately 
15  feet.  The  boring  denoted  by  "C"  was  explored  to  a  depth  of  90  feet  using 
5-inch  Shelby  tubes  sampling  at  5-foot  intervals.  The  "SW"  boring  was 
examined  utilizing  a  split  spoon  sampler  and  standard  penetration  tool.  No 
undisturbed  samples  were  obtained,  but  density,  moisture  content,  and  soil 
types  were  obtained  to  a  depth  of  90  feet.  The  "C"  hole  is  located  approxi¬ 
mately  in  the  center  of  the  test  area. 

Figure  XI-2  is  a  sum  ary  of  the  stress-strain  curves  obtained  from  the 
constrained  modulus  tests.  These  tests  were  performed  in  a  modified  Karol- 
Warner  Model  35^  consolidometer.  This  apparatus  provides  a  rise  time  to 
maximum  load  of  approximately  200  msec  and  essentially  zero  lateral  strains. 

The  results  show  composite  curves  for  all  the  samples  tested.  The  hearing 
capacity  of  the  surface  material  was  determined  from  triaxial  tests  to  be 
5,000  to  6,000  lbs /ft2. 

Investigations  of  various  sites  in  the  area  resulted  in  the  selection 
of  a  borrow  pit  in  the  area  of  the  experiment.  It  was  found  that  the  native 
material  could  be  placed  at  an  average  density  of  98  lbs/ft ^  by  conventional 
earth-moving  equipment  with  no  compactive  effort. 

it.  Predi  ctions 


Depth 


Seismic  velocity 


The  peak  vertical  stress  acting  at  any  depth  below  the  surface  of  the  loaded 
area  was  calculated  using  the  approximate  method  developed  by  Newmark  (reference  7). 
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This  method  utilizes  influence  charts  to  solve  the  classic  Boussinesq  spatial 
attenuation  problem  (i.e.,  a  distributed  static  load  applied  to  the  surface  of 
an  elastic  half-space).  Inherent  in  the  method  is  the  assumption  that  the 
stresses  propagate  at  an  infinite  velocity  and  thus  a  disturbance  at  the 
surface  immediately  influences  the  entire  body. 

The  idealized  pressure-time  curve  given  in  figure  1-13  was  utilized  as 
the  loading  function  with  the  smooth  curve  replaced  by  many  steps  of  uni  form 
pressure.  Each  step  contained  the  same  impulse  as  the  portion  of  the  curve 
that  it  replaced.  A  time  dependent  trial  and  error  solution  was  then  per¬ 
formed  with  the  load  positioned  at  various  points  over  the  96  x  150  foot  area 
until  the  maximum  stress  at  the  depth  of  interest  was  found.  This  process  was 
continued  until  the  maximum  stress  at  several  depths  was  found  and  then  the 
points  were  connected  with  a  smooth  curve. 

Figure  XI-3  presents  the  computed  theoretical  attenuation  of  maximum 
stress  with  depth  for  the  center  of  the  loaded  area.  The  results  for  the 
idealized  primacord  air-blast  input  and  a  i-MT  nuclear  air-blast  input  are 
shown  for  comparison. 

The  prediction  of  peak  vertical  acceleration  as  a  function  of  depth  is 
shown  in  figure  XJ-U.  Reference  7  was  used  as  the  basis  for  these  calcula¬ 
tions.  At  the  surface  the  peak  acceleration  was  computed  from  the  formula 


where 


Ap  «  peak  vertical  acceleration ,  g 

Pso  =  peak  overpressure  applied  at 
the  surface,  psi 

and 


C  =  seismic  velocity,  ft/sec 

At  depth,  the  peak  accelerations  were  computed  from  the  formula 
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** 


where 


peak  vertical  acceleration  at 
depth  of  interest,  g 


z  *  depth  of  interest,  ft 


a  “  attenuation  factor,  figure  XI -3 
.  z 

The  peak,  vertical  acceleration  at  a  depth  of  10  feet  was  connected  to  the 
surface  value  with  a  straight  line  to  interpolate  between  these  two  values. 


The  prediction  of  peak  vertical  velocity  as  a  function  of  depth  is  3hown 
in  figure  XI-5.  References  7  and  8  were  used  as  the  basis  for  these  calcula¬ 
tions.  Reference  7  gives  the  following  equation  for  computing  the  peak 
vertical  velocity 

Ifisao] 


V  =  U  ft/sec 
pz 


P 

so 

100 


where 


C  =  propagation  velocity,  ft/sec 

The  variations  in  results  are  obtained  by  using  different  values  for  C.  The 

change  in  seismic  velocity  at  58  feet,  as  noted  ir.  table  XI-1,  was  not  taken 

into  account .  The  lower  range  of  values  is  obtained  by  letting  C  equal  the 

seismic  velocity  of  1,700  ft/sec.  Hie  upper  range  of  values  is  obtained  by 

letting  C  equal  a  propagation  velocity  determined  from  the  constrained 

modulus  tests  (885  ft/sec).  It  should  be  noted  that  the  attenuation  with 

depth  is  determined  by  a  .  Reference  8  gives  the  following  equation  for 

z 

computing  the  peak  vertical  velocity  at  a  depth  of  5  feet. 

75Pso 

Y  -  — - —  ft/sec  +20  percent 


where 

V,.,  =  peak  vertical  velocity  at  a  depth 
J  of  5  ft ,  ft/sec 

s  **  specific  gravity 

*  3 A  of  the  seismic  velocity,  ft/sec 

This  calculation  contains  an  uncertainty  of  +20  and  this  entire  range 
has  been  indicated  on  figure  XI-5. 
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The  predicted  curve  was  drawn  approximately  midway  between  the  range  of 
vlaues  obtained  from  reference  7.  It  should  be  noted  that  the  predicted 
curve  also  agrees  quite  well  with  the  value  obtained  from  reference  8  for 
a  depth  of  5  feet. 

A  computer  program  was  available  to  calculate  both  relative  and  total 
displacements  for  various  depths.  This  program  is  patterned  after  the 
computational  scheme  presented  in  reference  7.  Tabulated  values  of  a 
and  P(t)  are  required  for  each  class  of  problems.  In  addition,  the  physical 
properties  of  the  medium  must  be  supplied.  Table  XI -2  contains  a  summary 
of  the  physical  properties  used  for  this  calculation. 

Table  XI -2 

PHYSICAL  PROPERTIES  FOR  DISPLACEMENT  CALCULATIONS 


Depth 

(ft) 

Density 

(lb/ft3) 

Seismic 
velocity 
( ft/sec) 

Wave  peak 
velocity 
( ft/sec) 

Recovery 

ratio 

(*) 

0-58 

no 

1,700 

885 

5 

58-268 

110 

2,550 

1,300 

100 

268-  ® 

no 

6,740 

3,400 

100 

The  wave  peak  velocity  is  the  velocity  at  which  the  peak  stress  propagates; 
this  value  should  be  determined  from  laboratory  tests.  For  the  deeper 
layers  the  wave  peak  velocity  was  assumed  to  be  approximately  one-half  of 
the  seismic  velocity.  The  recovery  ratio  is  a  measure  of  the  materials 
nonelasticity.  The  recovery  ratio  equals  the  recovered  strain  divided  by  the 
peak  strain  or 


R.R.“ 


seak  -  e  permanent  _  e  recovered 
e  peak  ~  e  peak 


(See  figure  XI -6.) 


Figure  XI-7  presents  the  prediction  of  total  displacements  as  a  function 
of  time  for  various  depths.  Figure  XI-8  presents  the  prediction  of  relative 
displacements  between  the  surface  and  various  depths  as  a  function  of  time. 
It  should  be  noted  that  the  majority  of  the  displacements  are  predicted  to 
occur  in  the  upper  75  feet  of  material. 
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The  angle  of  entry  of  the  stress  wave  into  the  ground,  a,  was  computed 
from  the  formula 


sin  a  *  y 


where 


U  -  air  shock  velocity,  U,80Q  ft/sec 


The  value  of  a  ranges  between  approximately  10  and  20  degrees  depending  on 
whether  C  equals  the  seismic  velocity  or  a  peak  stress  propagation  velocity , 
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XI-1.  Soil  Profile 
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Curves  from  Constrained  Modulus  Tests 
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Figure  XI-3.  Predicted  Attenuation  Factor 
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Figure  Xl-5.  Predicted  Peak  Vertical  Velocity 
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Figure  XI-8.  Predicted  Relative  Displacements  Between  Surface  and  Various  Depths 
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SECTION  XII 
PROCEDURE 


1.  Instrumentation 

The  instrumentation  used  in  this  experiment  was  selected  to  be  compatible 
with  the  recording  equipment  available  at  the  Air  Force  Weapons  Laboratory. 
Piezoelectric  devices  were  avoided  because  of  difficulties  with  long  lines 
and  with  moisture.  Further,  because  of  the  long  duration  of  the  loading 
pulse  in  the  experiment,  it  was  felt  that  piezoelectric  concepts  would  not  give 
sufficient  static  response  characteristics.  Wherever  possible,  balanced 
full -bridge  transducers  were  used.  DC  signal  conditioning  equipment  was 
preferred  over  carrier  systems . 

a.  Velocity  Gages 

The  primary  earth  free-field  measurement  in  this  experiment  was 
vertical  particle  velocity  at  various  depths.  The  velocity  gage  developed 
by  the  Sandia  Corporation  was  used  to  observe  this  motion.  This  gage  is  a 
3-Kc  carrier-driven  device  using  a  pendulum  blade  moving  through  a  viscous 
damping  fluid  to  provide  the  velocity  measurement.  The  gage  was  encased  in 
an  exterior  aluminum  shell  to  protect  it  from  the  soil  stress  environment . 
Figure  XII-1  shows  the  velocity  gage  and  protective  cover.  Calibration  was 
accomplished  in  the  laboratory  by  displacing  the  blade  pendulum  from  its 
zero  position  with  a  magnet  and  then  letting  the  1-g  gravitational  force 
return  the  pendulum  to  its  zero  position,  thus  producing  a  velocity  versus 
electrical  output  curve. 

b.  Accelerometers 

The  Statham  Model  A69TC  fluid  damped  accelerometer  was  used  to 
measure  vertical  accelerations  in  the  earth  free  field  at  various  depths. 

This  gage  was  selected  because  of  its  high  frequency  response.  The  accelero¬ 
meter  was  sealed  in  a  waterproof  aluminum  cup  to  protect  it  from  moisture  and 
soil  stress.  The  void  space  inside  the  aluminum  cup  was  filled  with  dry 
Ottawa  sand  to  improve  the  density  match  of  the  accelerometer  with  the  soil. 
Figure  XII-2  shows  the  accelerometer  and  protective  cover. 
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The  aceelerometers  were  calibrated  in  the  laboratory  on  a  centrifuge  prior 
to  packaging  for  placement  in  the  ground. 

c.  Soil  Stress  Gages 

A  newly  developed  soil  stress  gage  was  used  to  measure  the  vertical 
stress  in  the  soil  at  various  depths.  This  gage,  developed  at  the  Air 
Force  Shock  Tube  Facility,  utilizes  semiconductor  strain  gages  to  measure 
the  load  applied  to  a  small  aluminum  cylinder.  The  gage  was  calibrated  with 
static  air  pressure.  This  method  of  calibration  has  shown  excellent  agree¬ 
ment  with  results  when  the  gage  is  placed  in  a  dense  (112-113  lbs /ft3) 

Ottawa  sand  and  subjected  to  a  static  load.  No  dynamic  results  are  available. 
Figure  XII-3  shows  the  soil  stress  gage. 

d.  Long  Span  Displacement  Gage 

A  long  span  displacement  gage  was  used  to  measure  the  relative  dis¬ 
placement  between  the  surface  and  various  depths.  This  gage  was  developed 
and  has  been  previously  used  by  the  Sandia  Corporation.  The  device  consists 
of  a  long  wire  connected  to  concrete  deadmen  at  various  depths  in  a  drill 
hole.  The  wire  is  isolated  from  the  soil  by  flexible  metal  tubing.  The 
upper  end  of  the  wire  is  connected  to  a  spring-loaded  reel  whose  axle  is 
connected  to  a  rotary  potentiometer.  This  reel  is  attached  to  the  interior 
surface  of  a  cylindrical  concrete  structure  near  the  surface  of  the  ground. 

As  the  surface  is  depressed,  the  reel  moves  downward  and  the  rotation  of  the 
potentiometer  gives  an  electrical  signal  proportional  to  the  deflection. 

Figure  Xll-^  shows  the  long  span  displacement  gage  mounted  on  its  calibration 
device. 

e.  Time  of  Arrival  System 

A  modified  version  of  a  previously  developed  time-of-arrival  system 
was  used  to  monitor  the  movement  of  the  pressure  pulse  front  through  the 
earth.  The  system  consisted  of  100  onaii  directional  "ball  switches"  positioned 
in  a  three-dimensional  array  in  the  ground.  As  the  pressure  pulse  passes , 
the  switch  is  moved  and  gives  a  switch  closure.  These  switches  are  connected 
to  a  remote  console  where  the  switch  closure  triggers  a  multivibrator  which 
momentarily  lights  a  small  neon  light.  A  high-speed  camera  photographs 
the  light  hank  to  give  a  record  of  the  switch  closures  as  a  function  of  time. 
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f.  Slope-Indicator  Gage 

This  installation  was  used  to  measure  the  permanent  vertical  and 
horizontal  displacements  and  the  verical  transient  displacement  from  the 
ground  surface  to  a  depth  of  50  feet.  The  gage  consists  of  sections  of 
extruded  aluminum  casing,  2  feet  long  and  3  inches  in  diameter  with  instru¬ 
ment  grooves  90°  apart,  joined  together  with  12-inch  aluminum  couplings. 

The  couplings  permit  relative  vertical  movement  between  the  segments  and 
form  a  50-foot-long  composite  casing.  A  first-order  survey  was  performed 
pre-  and  post-shot  to  determine  the  absolute  notion  of  the  top  of  the  casing. 

The  change  in  inclination  of  the  casing  and  the  vertical  movement  of  the 
segments  relative  to  the  top  are  measured  with  an  inclinometer  or  slope 
indicator  instrument  and  depth  probe  lowered  down  the  inside  of  the  casing 
in  the  instrument  grooves.  These  measurements  were  made  pre-  and  post-shot; 
the  difference  between  the  two  sets  of  readings  represents  the  permanent 
horizontal  and  vertical  movement  of  the  casing. 

In  an  attempt  to  measure  vertical  transient  displacements,  a  50-foot 
section  of  continuous  blued  aluminum  tubing  2  1/2  inches  in  diameter  was 
inserted  inside  the  outer  casing  and  secured  to  the  bottom.  This  tubing 
was  installed  after  the  initial  set  of  inclinometer  and  depth  probe  measure¬ 
ments  had  been  made.  Spring-loaded  hinges  with  a  scratch  point  were  affixed 
to  the  inside  of  each  coupling  on  the  outside  casing  and  were  held  open  with 
a  continuous  length  of  piano  wire.  After  the  inside  blued  tubing  was  installed, 
the  piano  wire  was  withdrawn,  allowing  the  hinges  with  their  scribes  to  snap 
against  the  blued  tubing.  Figures  XII-5,  XII-6,  and  XII-7  illustrate  the 
slope  indicator  instrument  and  the  depth  probe. 

2.  Gage  Placement 

With  the  exception  of  three  near-surface  soil  stress  gages,  all  free- 
field  soil  transducers  were  placed  in  8-inch  diameter  drill  holes  at  various 
depths  ranging  from  1  foot  to  80  feet  depth  of  burial.  The  holes  were 
drilled  using  a  special  drilling  mud  which  formed  a  starch  liner  around  the 
interior  surface  of  the  holes.  This  liner  added  virtually  no  strength, 
but  prevented  water  from  penetrating  into  the  surrounding  soil  and  provided 
enough  stability  to  the  loose  granular  material  to  keep  the  holes  open 
until  the  instrumentation  was  placed. 
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Each  transducer  was  arranged  on  its  cable  so  as  to  be  properly  oriented 
when  it  was  lowered  into  the  drill  hole.  Short  sections  of  copper  tubing 
were  installed  over  the  soil  stress  gage  cables  to  ensure  vertical  align¬ 
ment.  In  addition,  a  light  was  lowered  into  the  hole  immediately  following 
the  gage  to  permit  visual  inspection  of  the  positioning  of  the  gage.  After 
the  placement  of  each  gage,  the  hole  was  filled  to  the  next  level  with  a 
dry  native  sand.  A  wire-screen  box  was  placed  over  the  top  of  the  hole  to 
provide  a  uniform  method  of  raining  the  sand  into  the  hole.  This  method 
gave  a  density  of  backfill  materiel  very  nearly  equal  to  that  of  the  surround¬ 
ing  soil.  The  three  ne sir-surface  soil  stress  gages  were  placed  in  hand-dug 
holes.  The  gage  holes  were  backfilled  with  a  native  sand. 

During  placement,  each  gage  was  connected  to  its  respective  signal 
conditioning  equipment.  All  gages  were  monitored  for  operation  and  balance 
to  ensure  their  proper  functioning. 

The  slope  indicator  casing  was  installed  in  a  9-inch-diameter  drill  hole 
and  backfilled  with  a  dry  native  sand.  The  sand  backfill  was  compacted  with 
a  pnuematic  vibrator  to  ensure  that  the  casing  would  closely  follow  the 
motion  of  the  surrounding  soil. 

3*  Gage  Locations 

Figure  XII-8  is  a  plan  view  of  the  test  area  showing  the  location  of 
all  of  the  earth  free-field  instrumentation.  It  will  be  noted  that  each 
position  containing  active  instrumentation  has  been  assigned  a  location 
number.  The  location  number,  gage  type  and  depth  of  burial  have  been  combined 
into  a  coding  system  which  is  illustrated  in  table  XII-1. 

Table  HI-1 

INSTRUMENTATION  CODING  SYSTEM 
Description 

Hole  19  -  long  Span  Displacement  -  75  ft  depth 
Hole  21  -  Soil  Stress  Gage  -  5  ft  depth 
Hole  25  -  Velocity  Gage  -  1  ft  depth 
Hole  12  -  Accelerometer  -  TO  ft  depth 
Hole  8  -  Soil  Stress  Gage  -  near  surface  (3  in.) 

Hole  6  -  Time  of  Arrival  -  20  ft  depth 


Code 

19D-75 

21S-5 

25V-1 

12A-70 

8S-g 

6T-2Q 
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Figure  XII-9  is  a  cross  section  of  the  typical  instrumentation  installa¬ 
tion.  Hole  21  contained  accelerometers,  velocity  gages,  and  soil  stress 
gages.  Holes  12  and  25  were  identical  to  hole  21  with  the  exception  that  all 
of  the  soil  stress  gage  locations  were  deleted.  The-  time-of-arrival  system 
was  normally  installed  in  either  a  25-  or  50-foot  deep  hole  as  illustrated 
in  figure  XEI-9.  Ihe  three  near  surface  (3  inches  deep)  soil  stress  gages 
were  installed  at  locations  8,  11,  and  16.  The  long  span  displacement  gage 


had  three  depths:  25 

feet,  50  feet,  and  75  feet. 

All  of  the  electronic 

instrumentation ,  with 

the  exception  of  the  time-of- 

•arrival  system,  responded 

to  only  the  vertical 

ccmpcnent  of  the  input. 

Table  XII-2  contains  a  complete  earth  free-field  instrumentation  listing. 

Table  XII-2 

INSTRUMENTATION  LISTING 

Hole  no. 

Gage  type 

Depths  ( ft ) 

2 

Time  of  Arrival 

5,10,15,20,25 

3 

Time  of  Arrival 

5,10,15,20  ,25 

6 

Time  of  Arrival 

1,5,10,15  ,20,25 

8 

Soil  Stress 

0.25 

9 

Time  of  Arrival 

5,10,15  ,20,25, 
30 ,35,^0,1*5,50 

11 

Soil  Stress 

0.25 

12 

Velocity 

1,20,5*0,60,80 

12 

Accelerometer 

10,30,50,70 

13 

Time  of  Arrival 

5,10,15,20,25, 

30,35,^0,1*5 

15 

Time  of  Arrival 

5,10,15 ,20,25, 
30,35,5*0,1*5 

16 

Soil  Stress 

0.25 

17 

Time  of  Arrival 

5,10,15,20,25 

19 

Long  Span  Displacement 

25,50,75 

19 

Accelerometer 

75 

20 

Time  of  Arrival 

5,10,15,20,25 

21 

Velocity 

1,20,1*0,60,80 

21 

Accelerometer 

10,30,50,70 
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Table  HI-2  (cont'd) 


Hole  no. 

Qa«e  type 

Depths  (ft) 

21 

Soil  Stress 

5,15,25,35, 

*♦5,55,65 

23 

Time  of  Arrival 

1,5,10,15,20,25, 
30, 35, U0 ,1*5,50 

25 

Velocity 

1,20,1*0 ,60,80 

25 

Accelerometer 

10,30,50,70 

26 

Time  of  Arrival 

5,10,15 ,20,25,30 
35,1*0,1*5 

29 

Time  of  Arrival 

5,10,15,20 ,25 

31 

Time  of  Arrival 

1,5,10,15 ,20,25 

30*35  * 

5,10,15 ,20,25 
1,5,10,15,20,25 


33 

35 


Time  of  Arrival 
Time  of  Arrival 


1 


I 
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Figure  XII-1.  Sandia  Velocity  Gage  and  Protective  Cover 
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Figure  XII-6.  Slope -Indicator  Instrument  and  Accessories 
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SECTION  XIII 
RESULTS 


1.  Soil  Stress  Gases 

Two  data  traces  were  recorded  from  the  10  soil  stress  gages  which  were 
installed  in  the  experiment.  All  the  soil  stress  gages  installed  in  Hole  21 
recorded  an  excessive  amount  of  noise  from  which  no  useful  data  could  be 
recognised  or  extracted.  The  noise  on  all  of  these  records  started  at 
approximately  13  msec,  which  agrees  with  the  first  signal  output  from  the 
top  velocity  gage  in  this  hole.  The  cause  of  this  excessive  noise  has  not 
been  positively  identified.  However,  it  is  suspected  that  the  noise  was 
generated  by  the  cables  being  squeezed  as  they  exited  from  the  protective 
conduit  at  the  drill  hole.  The  two  data  traces  were  obtained  from  gages 
installed  with  3  inches  of  earth  cover.  The  third  near-surface  gage  (llS-g) 
appeared  to  saturate  the  amplifiers  and  no  useful  data  were  obtained.  The 
measured  near-surface  vertical  soil  stress  records  for  gages  8S-g  and  l6S-g 
are  shown  in  figured  XIII-1  and  XIII-2  respectively. 

2.  Accelerometers 

Five  data  traces  were  recorded  from  the  13  accelerometers  which  were 
installed  to  measure  earth  free-field  motions.  The  remaining  eight  channels 
contained  an  excessive  amount  of  noise  from  which  no  useful  data  could  be 
recognized  or  extracted.  The  cause  of  this  excessive  noise  has  not  been 
positively  isolated;  however,  it  is  suspected  that  the  tape  recorders  contrib¬ 
uted  significantly  because  of  incorrect  gain  settings. 

Appendix  C  contains  the  time  history  records  for  all  of  the  accelerometers 
which  contained  data.  Figure  XIII-3  is  a  summary  of  the  recorded  peak  down¬ 
ward  vertical  accelerations  for  the  earth  free  field.  Acceleration  data 
from  the  interior  of  the  structural  model  have  also  been  included  on  this 
figure.  Figure  XIII-5  presents  the  times  of  first  arrival  of  the  accelera¬ 
tion  pulses  measured  in  Hole  12.  Figure  XIII-6  presents  the  time  of  first 
arrival  of  the  acceleration  pulse  measured  in  Hole  21. 

The  accele rat ion -time  records  were  integrated  to  obtain  additional 
transient  vertical  velocity  data  and  these  results  are  presented  in  appendix  C. 
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3.  Velocity  Gages 

Thirteen  data  traces  were  recorded  from  the  15  velocity  gages  which  were 
installed  in  the  experiment.  Gage  21V-1  was  disconnected  before  the  test 
because  of  a  faulty  signal  conditioning  unit.  Gage  21V-60  appeared  to 
saturate  the  amplifiers  and  no  useful  data  were  obtained. 

Appendix  C  contains  the  time  history  records  for  all  of  the  velocity 
gages  which  contained  data.  Figure  XIII-2+  is  a  summary  of  the  recorded 
peak  downward  vertical  velocities.  Figure  XIII-5  presents  the  times  of 
first  arrival  for  each  of  the  velocity  pulses  measured  in  Hole  12.  Figure 
XIII-6  presents  the  times  of  first  arrival  for  each  of  the  velocity  pulses 
measured  in  Hole  21.  Figure  XIII-7  presents  the  times  of  first  arrival  for 
each  of  the  velocity  pulses  measured  in  Hole  25.  A  summary  of  all  the  times 
of  arrival  of  the  peak  vertical  downward  velocity  are  presented  in  figure 

xiii-8. 

U.  Long-Span  Displacement  Gages 

The  long  span  displacement  gages  did  not  provide  any  useful  data.  The 
overpressure  caved  in  the  protective  cover  over  the  gage  housing  and  sub¬ 
jected  the  potentiometers  to  the  free-field  air  blast  environment.  This 
drove  all  of  the  displacement  records  entirely  off  scale  and  saturated  the 
amplifiers. 

The  velocity-time  records  were  integrated  to  obtain  an  indication  of 
the  transient  vertical  displacement;  these  results  sire  presented  in  appendix 
C.  Base  line  shifts  were  required  to  integrate  a  portion  of  the  velocity 
records  and  the  adjustments  that  were  made  are  indicated  by  dashed  lines  on 
the  figures  in  appendix  C.  A  summary  of  the  first  major  downward  peak 
vertical  displacements  obtained  from  the  integration  of  the  velocity  records 
is  presented  in  figure  XIII -9. 

5.  Slope-Indicator  Gage 

Figure  XIII-1Q  shows  a  plan  view  of  the  location  of  the  top  of  the  slope 
indicator  gage  both  pre-  and  post-shot  as  determined  from  survey  data.  These 
data  indicate  that  the  top  of  the  gage  moved  3.72  inches  downward,  1.20  inches 
to  the  east  and  O.lU  inch  to  the  south.  Figures  XI1I-11  and  XIII-12  illustrate 
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the  permanent  lateral  motion  of  the  gage  with  depth.  Table  XII I -1  gives  the 
permanent  vertical  displacement  with  depth. 

The  scratch  gages  installed  to  measure  the  transient  vertical  movement 
did  not  function  as  expected.  Problems  were  encountered  in  the  installation 
of  the  gage  and  the  magnitude  of  the  displacements  appear  to  be  lover  than 
expected.  These  factors  make  the  recognition  of  the  data  impossible  and, 
therefore ,  no  results  are  presented. 

Table  XIII-1 

PERMANENT  VERTICAL  DISPLACEMENT  WITH  DEPTH 
(SLOPE  INDICATOR) 

Displacement  Displacement  between 

Depth  surface  to  depth  successive  depths 

(ft)  (in.)  (in.) 


2.00 

1.32 

1.32 

1*.50 

1.80 

0.b8 

7.02 

2.16 

0.36 

9.52 

2.70 

0.2b 

12.0b 

2. bo 

0.00 

lb.  5b 

2.52 

0.12 

17.05 

2.70 

-0.12 

19.55 

2.52 

+0.12 

22.06 

2.52 

0.00 

2b.  57 

2.6b  (Max) 

+0.12 

27.08 

2.52 

-0.12 

29.58 

2.52 

0.00 

32.09 

2.6b  (Max) 

+0.12 

3b.  60 

2.52 

-0.12 

37.09 

1.32 

-0.12 

39.60 

2.28 

+0.96 

b2.09 

2. b0 

+0.12 

bb.62 

— 

— 

f 

l 
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6.  Tiiae-of-Arrival  System 


A  plan  viev  of  the  test  area  showing  the  location  of  the  tirae-of-arrival 
gages  and  the  shock  wave  at  the  surface  of  the  soil  is  presented  in  figure 
XIII-13*  As  indicated  in  the  figure,  the  fcime-of-arri  val  gages  were  installed 
along  three  lateral  lines.  The  results  obtained  from  the  time -of -arrival 
system  are  presented  in  figures  Xlll-lk,  XIII-15  ,  and  XIII-16  for  Planes  1,2, 
and  3,  respectively.  Each  plane  represents  an  east-vest  profile  of  the  test 
area.  The  time  of  arrival  at  each  gage  location  is  noted  on  the  figure, 
and  contour  lines  of  equal  times  have  constructed  to  indicate  the  location 
and  shape  of  the  stress  wave  in  the  earth. 
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Figure  XIII-1.  Near-Surface  Vertical  Soil  Stress 
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PEAK  VERTICAL  ACCELERATION.^'* 

Figure  XIII-3.  Peak  Vertical  Accelerations 
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Figure  XIII-4.  Peak  Vertical  Velocities 
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Figure  XIII-5.  Times  of  First  Arrival,  Hole  12 
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Figure  XIII-7.  Times  of  First  Arrival,  Hole  25 
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SECTION  XIV 
DISCUSSION 

1.  Soil  Stress 

As  indicated  in  figure  IX-1,  The  air-blast  input  had  a  different  shape 
and  considerably  faster  decay  to  one-half  peak  pressure  than  predicted.  The 
shape  of  the  air-blast  curve  influences  the  attenuation  of  stress  with  depth 
and,  therefore,  an  after-the-fact  calculation  of  the  stress  attenuation  was 
accomplished  utilizing  the  composite  air-blast  curve  given  in  figure  IX-1. 

The  results  of  this  calculation  sure  presented  in  figure  XIV-1  along  with  the 
predicted  curve.  Since  there  are  considerable  differences  in  the  two  curves, 
the  attenuation  of  other  earth  free-field  quantities  was  recomputed  after-the- 
fact  utilizing  the  new  attenuation  factors.  Complete  confidence  cannot  be 
placed  in  either  calculated  curve  without  supporting  stress  measurements , 
which  were  not  available  from  this  test. 

The  measured  vertical  stress  in  the  soil  near  the  ground  surface  (3  inches 
of  cover)  is  presented  in  figures  XIV-2  and  XIV-3.  These  records  indicate 
values  of  peak  overpressure  of  TOO  psi  and  570  psi  instead  of  the  312  psi 
obtained  from  a  normalization  of  all  the  air  pressure  measurements  (figure 
IX-l) .  The  composite  air  pressure  curve  has  been  plotted  on  each  of  the  soil 
stress  records  for  comparison.  The  overpressure  axis  has  been  scaled  such 
that  the  maximum  values  from  each  curve  coincide.  It  should  be  noted  that 
the  shape  of  the  two  curves  agrees  quite  well.  If  some  adjustments  were 
made  for  the  small  amount  of  attenuation  which  should  take  place  in  the  3-inch 
depth,  the  two  curves  would  agree  even  better.  Thus,  it  appears  that  the 
near  surface  stress  gages  reproduced  the  time  history  of  the  air-blast  input 
closely  while  magnifying  the  intensity  by  a  factor  of  approximately  2. 

These  gages  were  calibrated  with  static  air  pressure.  This  method  of  calibra¬ 
tion  has  been  correlated  with  results  in  a  dense  (112-113  lbs /ft3)  Ottawa 
sand,  subjected  only  to  static  loads.  The  gages  in  this  test,  however,  were 
placed  in  a  native  sand  with  a  density  of  approximately  100  lbs/ft3  and  were 
subjected  to  a  dynamic  load.  A  magnification  factor  of  2  is  not  unreasonable 
when  the  differences  in  density  and  loading  are  considered. 
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2,  Acceleration 

Figure  XIV-b  presents  all  of  the  measured  earth  free-field  peak  accelera- 
tions  along  with  after-the-fact  calculations  of  attenuation  with  depth,  which 
takes  into  account  the  composite  air  pressure  curve  (figure  IX-l)  and  the 
after-the-fact  attenuation  factors  given  in  figure  XIV-1.  Data  from  the 
interior  of  the  structural  model  have  also  been  included  on  this  figure.  The 
data  obtained  from  the  earth  free-field  accelerometers  were  of  extremely  poor 
quality  since  the  signals  were  masked  by  a  high-intensity  noise.  The  only 
record  for  which  this  was  not  true  was  19A-75  j  it  gave  a  peak,  acceleration 
which  was  much  higher  than  calculated.  All  of  the  acceleration  data  indicated 
pulses  of  extremely  short  duration.  When  these  data  were  integrated  to 
obtain  velocity-versus-time  histories,  the  results  indicate  that  none  of  the 
velocities  returned  to  zero  (appendix  C).  This  uncertainty  about  the  integrated 
results  also  raises  a  question  as  to  whether  the  peak  velocity  had  been 
attained  before  the  end  of  the  data. 

In  spite  of  all  the  uncertainties  ,  the  data  taken  as  a  whole  indicate 
that  the  shape  of  the  attenuation  of  peak  vertical  acceleration  curve  with 
depth  agrees  fairly  well  with  the  calculated  curve,  although  the  measured 
values  in  the  upper  20  feet  were  somewhat  lower.  A  much  larger  quantity 
of  reliable  data  would  be  required  before  the  prediction  technique  could  be 
radically  changed. 

3.  Velocity 

The  vertical  velocities  measured  at  the  near-surface  locations  had  a 
characteristic  shape  that  agreed  quite  well  with  the  air-blast  input,  i.e., 
a  short  rise-time  to  the  peak  followed  by  a  decay  to  zero  in  a  time  ranging 
from  approximately  100  to  200  msec  (see  figures  C-17  and  C-27  in  appendix  C) . 

The  rise  time  increased  monotonically  with  depth  as  indicated  in  table  XIV-1. 
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Table  XIV-1 
RISE  TIMES 


Hole 

Depth 

(ft) 

ta 

(msec) 

tp 

(msec) 

tp-ta 

(msec) 

12 

20 

3l* 

70 

36 

12 

i*o 

1*2 

85 

1*3 

12 

80 

1*8 

170 

112 

21 

1 

13 

ll* 

1 

21 

20 

31* 

85 

51 

21 

1*0 

1*2 

100 

58 

21 

60 

1*1* 

105 

61 

21 

80 

52 

120 

68 

25 

1 

..  12.2 

18 

5.8 

25 

20 

21*. 6 

50 

25.1* 

25 

1*0 

30 

65 

35 

25 

60 

31 

60 

29 

25 

80 

1*5 

110 

65 

NOTE:  ta  *  tine  of  first  arrival 
tp  *  time  of  peak  arrival 

In  addition,  the  measured  peak  vertical  velocity  generally  decreases  with 
depth  as  shown  in  figure  XIV- 5 .  This  figure  presents  all  the  measured  peak 
vertical  velocities  and  the  after-the-fact  calculated  attenuation  curve  as 
well  as  the  results  obtained  from  the  integration  of  acceleration  records. 
The  after-the-fact  calculations  were  based  on  the  formula  from  reference  7, 


which  is  given  in  section  XI 


V 

pz 


A  peak  over- 


( 

{ 


l 

i 


i 


t 

pressure  of  312  psi  was  used  for  P  ,  C  was  taken  as  3/4  of  the  seismic  l 

SO  I 

values  given  in  table  XI -1,  and  o  was  taken  from  the  after-the-fact  calcula-? 

z  ! 

tions  given  in  figure  XIV-1.  Although  there  is  considerable  scatter,  ‘ 

particularly  when  the  velocities  obtained  by  integration  are  included,  the  { 
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data  appear  to  fit  the  calculations  reasonably  well  with  the  exception  of  the 
near  surface  values.  Reference  8  does  not  predict  velocity  values  above  a 
depth  of  5  feet  because  of  the  erratic  results  observed  from  nuclear  field 
tests.  The  technique  for  predicting  accelerations  given  in  reference  7  takes 
this  into  account;  near  surface  values  and  surface  values  are  computed  by 
two  different  equations  which  give  widely  different  results.  It  is  also 
possible  that  the  velocities  are  much  higher  at  the  surface  than  at  small 
depths  of  burial.  This  could  be  caused  by  low  initial  densities  of  the 
surface  material  or  by  a  scouring  effect  immediately  ahead  of  the  air  blast. 

The  data  indicate  that  the  velocity  at  a  depth  of  20  feet  is  much  higher  in 
Hole  t.l  (9*9  ft/sec)  than  in  Holes  12  (4.7  ft/sec)  and  25  (5-8  ft/sec).  A 
first  reaction  would  be  to  assume  that  this  efeect  was  caused  by  differences 
in  the  -Jr  pressure  input  at  these  locations.  A  review  of  the  pressure-time 
records  ^iven  in  appendix  B  reveals  that  there  was  no  record  obtained  from 
the  gage  near  Hole  21.  Both  of  the  air  pressure  gages  (10  and  24)  near  Holes 
12  and  indicate  that  the  pressure  records  required  base  line  shifts  to 
correct  for  negative  pressures  at  relatively  early  times.  In  addition,  the 
air  pressure  record  from  Position  24  was  very  spiked  in  appearance.  However, 
there  were  not  sufficient  air  pressure  measurements  taken  to  positively 
determine  the  variation  in  air  pressure  input  at  various  locations  within  the 
loaded  area.  In  addition,  the  near  surface  velocities  had  fair  agreement, 
indicating  similar  air-blast  inputs. 

Another  possible  cause  of  these  differences  might  be  the  side  rarefactions 
as  illustrated  in  figure  1-5.  However,  this  does  not  appear  to  be  the  case. 
Figure  XIV-6  illustrates  a  gross  simplification  of  the  real  problem,  i.e.,  a 
static  load  applied  by  a  foundation  to  a  soil  which  is  assumed  to  have  elastic 
properties.  In  this  case,  the  rarefactions  have  all  interacted  and  static 
equilibrium  has  been  reached.  The  unit  width  can  be  used  to  represent 
approximately  100  feet.  Then  the  centerline  represents  Hole  21  and  the 
one-quarter  point  represents  Holes  12  and  25.  At  a  scaled  depth  of  20  feet 
the  stress  intensity  is  90  percent  and  94  percent  respectively.  Thus  the 
static  representation  of  the  rarefaction  phenomenon  does  not  explain  the 
6 0  percent  or  more  differences  in  the  measured  velocities.  Another  argument 
can  be  used  to  negate  the  importance  of  the  rarefaction  waves.  Figure  XIV-7 
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illustrates  the  geometry  and  travel  times  involved.  In  this  figure,  t  *  time 

of  first  arrival  for  the  signal  generated  at  the  surface,  t  *  time  of  peak 

arrival  for  the  signal  generated  at  the  surface,  and  t  *  time  of  first 

r 

arrival  for  the  rarefaction  wave  generated  at  the  boundary  of  the  loaded 
area.  The  fastest  first-arrival  velocity  was  utilized  to  compute  the  arrival 
times  for  the  rarefaction  waves.  It  is  obvious  from  the  figure  that  the 
rarefaction  waves  arrive  at  all  holes  after  the  first  arrival  from  the  surface 
but  before  the  peak  arrival  from  the  surface.  The  arrival  of  a  rarefaction 
wave  cannot  be  observed  on  the  velocity  records  given  in  appendix  C  between 
the  first  arrival  and  the  peak  arrival  for  the  gages  in  question. 


I 


Another  explanation  for  the  differences  in  velocity  measurements  at  the 
20-foot  depth  is  the  difference  in  the  in-situ  material  properties.  The 
stratigraphy  in  Hole  21  consists  of  10  feet  of  fractured  hard  clay  followed 
by  a  thick  layer  of  sand  to  a  depth  of  35  feet.  The  stratigraphy  in  Holes 
12  and  25  consists  of  a  foot  or  more  of  a  very  loose  silt,  a  4  to  5-foot 
layer  of  hard  fractured  clay  and  a  thick  layer  of  sand  (see  figure  XI-l), 

It  is  also  obvious  from  the  blow  counts  on  this  figure  that  there  were 
considerable  variations  in  the  density  of  the  material  in  the  sand  layer. 

Thus  the  measured  velocities  at  a  depth  of  20  feet  may  have  been  different 
because  of  differences  in  attenuation  which  took  place  above  this  depth  and/or 
because  of  differences  in  the  density  of  the  material  in  which  the  gages 
were  located. 


The  time -of- arrival  data  from  the  velocity  records  can  be  used  to  compute 
first-arrival  propagation  velocities  as  indicated  in  figures  XIV-8,  XIV-9, 
and  XIV-10.  The  propagation  velocities  computed  from  the  first  arrivals 
agree  with  the  uphole  seiBmic  data  (table  XI-l)  with  two  exceptions.  The 
propagation  velocity  in  the  upper  20  feet  of  material  in  Hole  21  is  approxi¬ 
mately  50  percent  of  the  seismic  value.  In  addition,  both  holes  21  and  25  : 

i 

indicate  an  exceedingly  fast  layer  between  20  and  1*0  feet  with  propagation 
velocities  in  the  range  of  10,000  to  20,000  ft/sec.  Referring  to  the  Soil  s 

Profile  (figure  XI-l),  the  high  velocity  region  appears  to  correspond  to  [ 

the  very  dense  sand  and  gravel  layer.  It  is  very  doubtful  that  &  dense  s 

sand  and  gravel  layer  would  be  capable  of  providing  velocities  of  this  j 

I 

magnitude.  Hie  hard  caliche  in  the  lower  layer  should  have  a  much  faster  j 
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velocity  than  the  sand  and  gravel,  but  the  results  again  show  good  agreement 
with  the  seismic  data. 

The  time-of-arrival  data  from  the  velocity  records  can  also  be  used  to 
compute  peak  arrival  propagation  velocities  as  indicated  on  figure  XIV-11. 
These  propagation  velocities  indicate  that  the  peak  velocity  travels  at  one- 
half  the  front  velocity  or  less.  This  agrees  quite  well  with  empirical 
rules  of  thumb  ranging  from  1/2  to  3 A  of  the  front  velocity. 

Displacement 

No  measured  data  were  obtained  on  transient  displacement  since  both  the 
long  span  displacement  gage  and  the  transient  measuring  portion  of  the  slope 
indicator  failed  to  perform  satisfactorily.  The  peak  vertical  displacements, 
derived  from  the  integration  of  the  velocity  records ,  are  compared  on  figure 
XIV-12  with  after-the-fact  calculations.  The  agreement  is  very  poor  and  the 
data  are  very  inconsistent.  A  detailed  study  of  the  displacement  time 
records  in  appendix  C  will  reveal  that  the  majority  of  these  data  do  not 
have  the  expected  clasic  shape  as  illustrated  in  figure  Xl-7.  The  peak 
displacements  from  records  which  have  the  more  classic  shape  also  tend  to 
fit  the  calculated  curve  on  figure  XIV-12  the  best.  These  discrepancies 
may  be  caused  by  inaccuracies  in  the  velocity-time  records  or  by  a  poor 
selection  of  base-line  shifts  which  were  utilized  for  the  integration. 

Permanent  displacement  data  are  available  from  three  sources :  inte¬ 
grated  velocity-time  records,  slope  indicator  measurements ,  and  pre-  and 
post-shot  surveys.  The  permanent  vertical  displacements  obtained  from  the 
integration  of  the  velocity  records  were  so  inconsistent  that  the  results 
are  not  presented.  The  permanent  displacement  of  the  top  of  the  slope 
indicator  installation  was  3.72  inches  according  to  the  results  obtained 
from  the  pre-  and  post-shot  surveys.  However,  the  depth  probe  readings 
indicate  a  value  of  2.6k  inches.  These  differences  in  results  have  not 
been  reconciled.  The  survey  also  revealed  that  all  of  the  near-surface 
monuments,  such  as  the  air  pressure  gage  mounts,  had  essentially  no  permanent 
displacement.  This  can  only  be  accounted  for  by  elastic  rebound  near  the 
free  surface  which  would  not  affect  deep  structures  to  the  same  extent.  The 
permanent  lateral  movement  in  the  east-west  direction  was  identical  for  the 
slope  indicator  and  the  survey  data.  The  results  indicate  the  surface 
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moved  1.2  inches  to  the  cast..  The  permanent  lateral  movement  in  the  north- 
south  direction  varied  between  the  slope  indicator  and  the  survey  data.  This  j 

f 

difference  appears  to  be  due  to  an  error  in  a  slope  indicator  reading  at  a  j 

j 

depth  of  1*6.5  feet  as  shown  in  figure  XIII-12.  A  permanent  lateral  movement  l 

l 

to  the  east  was  also  observed  in  the  structural  model.  A  movement  in  the  f 
direction  opposite  to  the  airblast  travel  was  not  expected  since  the  stress  \ 
wave  induced  in  the  ground  will  have  a  horizontal,  component  whose  direction 
is  to  the  west.  A  surveying  error  can  be  immediately  ruled  out  because  of 
the  excellent  agreement  between  the  two  independent  measuring  techniques  as 
illustrated  in  figure  XI II -11 .  This  motion  can,  however,  be  explained  by  an 

elastic  rebound  which  overshot  the  initial  position  during  the  unloading.  I 

* 

5.  Horizontal  Motions  and  Stresses  | 

? 

The  time-of-arrival  system  indicated  stress  wave  angles  of  entry  (see  | 
figure  1-7)  ranging  from  10°  to  1*5°.  However,  the  majority  of  the  angles  | 

ranged  from  15°  to  25°,  and  20°  will  be  taken  as  a  representative  average.  s 
Combining  the  equations  given  in  section  I  for  the  two  components  of  motion  f 
gives  | 


Horizontal  component  m  sin  a  _ 

Vertical  component  cos  a  8X1  a 

and  for  20°,  tan  a  *  0.361*.  Thus  the  Initial  horizontal  component  of  both 
motions  and  stresses  behind  the  wave  front,  is  equal  to  approximately  36 
percent  of  the  vertical  component. 
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Figure  XIV-1.  Attenuation  Factors 
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Figure  XI V-2 .  Wear-Surface  Vertical  Soil  Stress 


Figure  XIV-3.  Near-Surface  Vertical  Soil  Stress 
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Figure  XIV-1*.  Peak  Vertical  Accelerations 
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Figure  XIV-7.  Influence  of  Rarefaction  Waves 
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Figure  XIV-8.  Propagation  Velocities ,  Hole  12 
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Figure  XIV-9.  Propagation  Velocities,  Hole  21 
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Figure  XIV- 12.  Peak  Vertical  Displacements 
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SECTION  XV 


CONCLUSIONS  AND'  RECOMMENDATIONS 

The  peak  vertical  stresses  produced  "by  this  experiment  are  presented  in 
figure  XV-1.  Since  no  stress  data  were  obtained  at  depth,  this  curve  was 
prepared  utilizing  the  after-the-fact  attenuation  factors  given  in  figure 
XIV-1.  The  peak  vertical  accelerations  and  velocities  produced  by  this 
experiment  are  presented  respectively  in  figures  XV-2  and  XV-3.  The  vertical 
velocities  are  the  most  reliable  of  all  the  transient  earth  free-field 
phenomena  measured.  These  data  reflect  the  nonhomogeneous  nature  of  the 
soils  which  underlie  the  experimented  area.  The  peak  vertical  displacements 
produced  by  this  experiment  were  the  most  erratic  of  all  the  free-field 
quantities,  and  a  great  deal  of  reliance  cannot  be  placed  in  the  data.  It  is, 
therefore,  recommended  that  the  calculated  curve  presented  in  figure  XIV- 12 
be  utilized  for  defining  the  variation  in  peak  vertical  displacement  with 
depth . 

The  ground  surface  experienced  a  permanent  horizontal  displacement  of 
approximately  1  inch  in  a  direction  opposite  to  air-blast  propagation,  and 
a  negligible  amount  of  permanent  vertical  displacement.  However,  near-surface 
structures  experienced  permanent  vertical  displacements  ranging  from  2  to  It 
inches.  The  magnitude  and  direction  of  the  permanent  displacements  were 
primarily  governed  by  elastic  rebound  phenomena. 

The  selection  of  gages  for  this  experiment  as  regards  type,  number,  and 
location  should  have  been  adequate  to  meet  the  project  objectives  if  a  high 
degree  of  data  recovery  had  been  realized.  Mechanical  and  structural  changes 
must  be  made  to  both  the  long-span  displacement  gage  and  the  slope-indicator 
gage  to  successfully  measure  transient  displacements.  Electronic  and 
recording  problems  must  be  corrected  if  the  earth  free-field  accelerations 
and  stresses  are  to  be  completely  described. 

The  mass  of  material  realistically  loaded  by  this  simulation  technique 
is  a  function  of  both  the  size  of  the  loaded  area  and  the  shape  of  the  air- 
blast  input.  As  a  rule  of  thumb,  the  width  of  the  loaded  surface  should  be 
twice  the  depth  of  interest,  and  at  this  depth,  the  width  of  realistic 
simulation  is  limited  to  one-half  the  width  of  the  loaded  surface.  This  is 
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illustrated  in  figure  XV-4.  In  addition,  the  actual  air-blast  input  must 
siffiilate  the  overpressure-time  history  from  the  nuclear  yield  of  interest 
as  described  in  section  I.  These  considerations  should  produce  a  simulated 
8 tress  environment  which  does  not  deviate  more  than  20  percent  from  the 
desired  nuclear  environment. 
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Figure  XV- 2.  Peak  Vertical  Accelerations 
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Figure  XV-3.  Peak  Vertical  Velocities 
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APPENDIX  A 

INSTRUMENTATION  EQUIPMENT  LIST 


1.  Tape  Recorder 


a.  Make :  Ampex 
Model :  CP-100 

Ampex  Corporation  Instrumentation  Products 
934  Charter  Street,  Redwood  City,  Calif 

List  of  Characteristics: 

(1)  Tape  Transport 

(a)  Tape  Speeds:  60,  30,  15,  T  1/2,  3  3/4,  17/8  ips  standard 

(b)  Tape  Speed  Deviation:  +G.35  percent  maximum 

(c)  Start  Time:  5  seconds  or  less  at  60  ips  or  30  ips 

3  seconds  or  less  at  15  ips  or  7  1/2  ips 
1.5  seconds  or  less  at  3  3/4  ips  or  1  7/8  ips 

(d)  Stop  Time:  Maximum  of  2.0  seconds 

(2)  Direct  Record/Reproduce  System 

(a)  Frequency  Response  and  RMS  Sign&l-to-Noise  Ratio: 


RMS  Signal  to  RMS  Noise 
Bypass  Un filtered 

Tape  speed  Bandwidth  filtered  wideband 

( ips )  ( cps )  ( db )  ( db ) 


60 

300 

to 

250,000  +  3 

db 

30 

25 

30 

200 

to 

125,000  +  3 

db 

30 

24 

15 

200 

to 

60,000  +_  3 

db 

30 

23 

7  1/2 

200 

to 

25,000  +  3 

db 

28 

21 

3  3/4 

200 

to 

12,500  +  3 

db 

28 

19 

1  7/8 

200 

to 

6,250  +  3 

db 

28 

19 

(b)  Input  Consumption :  Approximately  1375  recording  level; 
operable  from  0.7  to  10.0  volts  rms . 

(w)  Input  Impedance:  Minimum  18,000  ohms  resistive,  in 
parallel  with  275  mi cromicro farads ,  unbalanced  to  ground. 

(d)  Output  Level:  1.0  volt  rms  nominal  across  a  10,000  ohms 
or  greater  impedance. 
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(e )  Output  Impedance:  Less  than  100  ohms . 

(f)  Control  Track  Generator:  Subcarrier  frequency  17  kc  or 
18.24  kc.  Modulating  frequency  60  eps  +0.02  percent.  Modulation  50  percent 
+5  percent. 

(3)  Power  Requirements 

(a)  Voltage:  105  to  125  volts,  single  phase,  48  to  62  cps  or 
380  to  420  cps  AC;  or  210  to  250  volts,  single  phase,  48  to  62  cps  or  380 
to  420  cps  AC;  or  26  to  30  volts  DC,  ripple  2  volts  pesk-tc-peak  maximum. 

(b)  Power  Consumption:  Approximately  375  watts  for  a  14 

track  system. 

(4)  Environment 

(a)  Temperature ;  Operating  +40°F  to  +125°F.  Storage  non¬ 
operating  -20°F  to  l60°F. 

(b)  Altitude:  Operating  10,000  feet;  nonoperating  20,000  feet. 

(c)  Relative  Humidity:  Up  to  95  percent  without  condensation 

space . 

b.  Make:  Consolidated  Electrodynamics 
Model:  PR-3300 

Consolidated  Electrodynamics  Corporation 

360  Sierra  Madre  Villa,  Pasadena,  Calif 

List  of  Characteristics: 

(1)  Transport 

(a)  Tape  Speeds:  Six  standard  speeds  in  3  pairs  (60,  30  ips; 
15,  7  1/2  ips;  and  3  3/4,  1  7/8  ips) 

(b)  Tape  Speed  Accuracy:  +0.25  percent 

(c)  Start  Time:  Less  than  3  seconds  at  a  speed  of  60  ips  and 
2  seconds  at  lower  speeds, 

(d)  Stop  Time:  Less  than  2  seconds  at  60  ips  and  1  second  at 


all  other  capstan  controlled  speeds. 

(2)  Direct  Record  and  Reproduce  System 

(a)  Frequency  Response  and  Signal -to-Noise  Ratio: 


Tape  Speed 
(ips) 

60 

30 

15 

7  1/2 

3  3/4 

1  7/8 

Bandwidth 

300  cps 

200  cps 

100  cps 

100  cps 

100  cps 

100  cps 

measured 

to 

to 

to 

to 

to 

to 

100  kc 

50  kc 

25  kc 

12.5  kc 

6.25  kc 

3.125  kc 
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Frequency 
response 
ref  to  1  kc 

♦3 

db  +3  db 

+3  db 

+3  db 

t?  db 

+3  db 

S-N-R(RNE) 
300  cps  to 
band  edge 

-35 

db  -35  db 

-35  db 

-35  db 

-35  db 

-35  db 

(b) 

Input  Level: 

1  volt  rms 

nominal 

(0  db)  to  produce  normal 

recording  level. 

(c)  Input  Sensitivity:  0.25  to  10  volts  nas. 

(d)  Input  Impedance:  20k  ohms ,  unbalanced  to  ground. 

{e)  Output  Level:  1  volt  rms  nominal  {0  db)  across  a  600-ohm 
load  impedance. 

( f )  Output  Impedance:  Less  than  100  ohms  in  series  with  25  mf 
unbalanced  to  ground. 

(g)  Distortion:  1  percent  3rd  harmonic  distortion  of  a  one  kc 
signal  and  less  than  0.6  percent  intermodulation  distortion  when  recording 
and  reproducing  at  normal  record  level  (0  db)  at  60  ips. 

2 *  Timing  Control  System 

a.  Make:  Hyperion 

Model:  Hi-ll*0-12-S 

Hyperion  Industries ,  Inc. 

Watertown,  Mass 

List  of  Characteristics 

(1)  WWV  Synchronize:  Automatically  synchronizes  time  code  generator 
with  WWV  with  accuracy  of  1  msec. 

(2)  Decimal  time -of- day  display;  continuous  decimal  in-line  display 
synchronized  with  output  time  codes.  Time  of  day  displayed  in  hours, 
minutes,  seconds,  tenths  of  seconds. 

(3)  Time  Base:  Internal  1  me  crystal  controlled  oscillator  with 

o 

stability  of  five  parts  in  10  per  day. 

(k)  Serial  Time  Code  Outputs: 

(a)  I  RIG- A 

(b)  IRIG-B 

(5)  Serial  Modulated  Time  Code  Output: 

0  to  5  volts  P-P  into  1  k 

IRIG  Modulated  into  100  kc  carrier 

IRIG  Modulated  into  1  kc  carrier 
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(6)  DC  Level  Shift  Output:  DC  level  shift  code  output  is  0  to  -6 
volts  P-P  into  1  k  ohm.  Rise  time  2  microseconds. 

(7)  Power  Requirements:  105  to  125  V,  60  to  MiO  cps. 

(8)  Operating  Environment:  -25°C  to  +55°C  in  95  percent  humidity. 

3-  Signal  Conditioning  Equipment 

a.  Make:  Consolidated  Electrodynamics  Corporation 

Model:  CEC  System  D  (Type  1-113B  Carrier  Amplifier) 

Consolidated  Electrodynamics  Corp. ,  Data  Recorders  Division 
360  Sierra  Madre  Villa,  Pasadena,  Calif 

List  of  Characteristics: 

(1)  Full  Scale  Output:  +J>  ma  into  a  26-ohm  load. 

(2)  Sensitivity:  Input  signal  for  maximum  output,  1  mv  with  no 
attenuation,  1  volt  maximum  with  full  attenuation. 

(3)  Input  Impedance:  Approximately  1,800  ohms. 

(4)  Input  Attenuator:  1  to  1,000  in  20  steps. 

(5)  Bridge  Balance  System:  Will  accommodate  four  external  bridge 
anas  conposed  of  wire  strain  gages  or  other  resistance  elements  or  two-arm 
variable -reluctance  transducers  suitable  for  use  at  3  kc. 

(6)  Linearity:  Output  current  proportional  to  input  voltage  within 
2  percent  of  maximum  output. 

(7)  Frequency  Response:  Galvanometer  trace  amplitude  constant 
(+2  percent)  for  modulating  frequencies  from  0  to  600  pcs. 

b.  Make:  Allegany 

Model:  Sensor  Analogue  Module  ( SAM-1 ) 

Allegany  Instrument  Company,  Division  of  Textron  Electronics,  Inc. 
1091  Wills  Mountain,  Cumberland,  Maryland 

List  of  Characteristics: 

(1)  Gain  Range:  100  to  2,000,ccontinuous  with  vernier 

(2)  Amplifier  Balance  Control  Range:  +100  microvolts  referred  to 
the  input. 

(3)  Balance  Range:  +10  percent  of  scale  standards. 

( 4 )  Linearity:  0.02  percent  of  full  scale  at  DC. 

(5)  Frequency  Response:  DC  to  40  kc.  Down  3  db  at  1»0  kc  with 
vernier  fully  counterclockwise. 

(6)  Noise:  7  microvolts  rms ,  DC  to  40  kc.  Two  microvolts  peak-to- 
peak,  DC  to  3  cps. 
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(7)  Drift:  0.3  ndcrovolt/°F  referred  to  input. 

(8)  Output  Impedance:  Less  than  0.1  ohm  at  DC. 

(9)  Output  Voltage:  10  volts  DC  or  peak  AC  to  20  kc. 

(10)  Output  Current:  100  ma  DC. 

(11)  Common  Mode  Rejection:  90  db  at  60  cps. 

4.  Air  Pressure  Gage 


a.  Make :  Norwood 

Model:  111  Bonded  Strain  Gage  Pressure  Transducer 

Advanced  Technology  Laboratories,  A  Division  of  American -Standard 
360  Whisman  Road,  Mountain  View,  Calif 

List  of  Characteristics: 

(1)  Nonlinearity:  Better  than  0.5  percent  of  F.S.  by  terminal 
method  or  0.25  percent  by  best  straight  line  through  zero  method. 

(2)  Hysteresis:  Better  than  0.5  percent  of  F.S. 

(3)  Repeatability:  Better  than  0.1  percent  of  F.S. 

(b)  Resolution:  Infinite. 

(5)  Acceleration  Effect:  Less  than  0.1  percent  F.S.  per  &  in  all 

planes. 


(6)  Vibration  Effect:  Insensitive  from  50  to  2,000  cps  to  100  & 
in  3  planes. 

(7)  Zero  Pressure  Output:  Less  than  +2  percent  of  F.S.  at  rated 
excitation. 


(8)  Resonant  Frequency:  Approximately  1*5,000  cps. 

(9)  Frequency  Response:  Flat  within  +  db  from  0  to  20,000  cps. 

(10)  Pressure  Limit:  150  percent  of  F.S.  for  static  pressure  (125 
percent  for  100  psi  range).  Full  scale  for  dynamic  pressures. 

(11)  Recommended  Excitation:  10  volts  DC  or  AC,  17  volts  maximum. 

(12)  Electrical  Output:  3.5  mv/v  +20  percent  for  -3l»  bridge,  3.0 
mv/v  +2  percent  for  -35  and  -36  bridge. 

b.  Make:  Scheavitz-Bytrex 
Model :  HFG-2000 

Bytrex  Corporation,  Kulite-Bytrex  Corporation 
50  Hunt  ST.i<rJt,  Newton,  Mass 

List  of  Characteristics: 

(l)  Output  Signal:  100  mv  F.S.  minimum. 
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(2)  Output  Impedance:  1,300  ohms  at  80°F,  increases  approximately 
0.015  percent/°F. 

(3)  Input  Impedance:  5C0  ohnss  at  80°F,  increases  approximately 

\ 

0.06  percent/°F. 

(4)  Excitation:  20  volts  DC  or  AC  rms. 

(5)  Natural  Frequency:  100  kc. 

(6)  Shock:  Will  Withstand  more  than  2,500  &  on  all  axes. 

(?)  Acceleration  Sensitivity:  0.002  percent  per  £. 

(8)  Operating  Temperature  Range:  -65  to  +300°F. 

(9)  Nonlinearity  and  Hysteresis:  Less  than  I  percent,  combined. 

(10)  Maximum  Allowable  Pressure:  150  percent  of  range  without  damage, 
200  percent  of  range  without  rapture . 

(11)  Measuring  Range:  0  to  2,000. 

5.  Accelerometers 

a.  Make :  Statham 

Model:  A69TC-350 

Statham  Instruments,  Inc. 

12U01  West  Olympic  Blvd,  Los  Angeles ,  Calif 

List  of  Characteristics: 

(1)  Acceleration  Range:  +100  £  to  jj>00 

(2)  Bridge  Resistance:  350  ohms. 

(3)  Excitation:  5  volts  DC  or  AC  rms. 

(H)  Full  Scale  Output:  Approximately  +20  mv. 

(5)  Operating  Temperature  Range:  -U0°F  to  +200°F. 

(6)  Direction  of  Sensitivity:  Perpendicular  to  base. 

(7)  Overload:  Three  times  rated  range. 

(8)  Transverse  Acceleration  Response:  Less  than  0.02  g/&  for 
transverse  accelerations  up  to  rated  range. 

(9)  Nonlinearity  and  Hysteresis:  Less  than  +1  percent  of  F.S, 

output . , 

(10)  Natural  Frequency:  Approximately  2,000  cps. 

6.  Soil  Pressure  Gages 
a.  Make :  Lynch 

Model:  Developed  and  fabricated  in  Air  Force  Weapons  Laboratory 
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List  of  Characteristics: 

(1)  Frequency  Response:  Greater  than  100  kc. 

(2)  Range:  0  to  greater  than  2,000  psi. 

7.  Velocity  Gage 

a.  Make:  Sandia  Corporation 

Model:  DX-B,  variable  reluctance  type 
List  of  Characteristics: 

(1)  Frequency  Response:  0.02  to  400  cps „ 

(2)  Maximum  Range:  100  fps. 

(3)  Excitation:  4  volts  rms. 

8-  Tlme-of-Arrival  System 

100-channel  digital  display  actuated  by  omnidirectional  inertia  switches, 
read  out  on  l6-mm  film  run  at  1,000  fps.  Supplies  a  three-dimensional 
analog  of  energy  transmission  in  soil. 

9-  Long-Span  Strain  Gage 

a.  Make:  Developed  and  fabricated  in  Air  Force  Weapons  Laboratory 
Model:  Helical  potentiometer  type 
List  of  Characteristics: 

(1}  Frequency  Response:  12  cps. 

(2)  Rise  Time:  21  msec. 

(3)  Range:  Up  to  100-ft  length. 
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Figure  B-l.  Air  Pressure  versus  Time 


"TF-  ■ 

- 


AFWL  TR-65-26,  Vol  I 


msec 


Figure  B-4.  Air  Pressure  versus  Time 
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Air  Pressure  versus 
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Figure  B-6.  Air  Pressure  versus 


Figure  B-7.  Air  Pressure  versus  Time 


Figure  B-8.  Air  Pressure  versus  Time 
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Figure  C-l.  Vertical  Acceleration  versus  Time 
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Figure  C-2.  Vertical  Velocity  versus  Time 
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Figure  C-5*  Vertical  Acceleration  versus  Time 
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Figure  C-6.  Vertical  Velocity  versus  Time 
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Figure  C-9.  Vertical  Acceleration  versus  Time 
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Figure  C-10.  Vertical  Velocity  versus  Time 
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Figure  C-13.  Vertical  Velocity  versus  Time 
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Figure  C-lU.  Vertical  Displacement  versus  Time 
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rtical  Displacement  versus  Time 
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Figure  19.  Vertical  Velocity  versus  Time 
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Figure  C-2.1 .  Vertical  Velocity  versus  Time 
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Figure  C-22.  Vertical  Displacement  versus  Time 
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Figure  26.  Vertical  Displacement  versus  Time 
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Figure  C-29.  Vertical  Velocity  versus  Time 
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Figure  C-30.  Vertical  Displacement  versus  Time 
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Figure  C-32.  Vertical  Displacement  versus  Time 
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Figure  C-33.  Vertical  Velocity  versus  Time 
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DISTRIBUTION 

HEADQUARTERS  USAF 
Hq  USAF,  Wash,  DC  20330 
( AFCOA) 

(AFOCE) 

(AMIN) 

{ AFRDPF ,  M&j  Dunn) 

USAF  Dep ,  The  Inspector  General  (AFIDI),  Norton  AFB,  Calif  92409 
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